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Analysis  of  material  damping  and  optimization  of  both 
material  damping  and  specific  stiffness  of  laminated, 
continuous  or  discontinuous  fiber  reinforced  polymer 
matrix  is  the  major  objective  of  this  study. 

Two  different  approaches,  laminated  plate  theory 
approach  and  energy  approach,  are  used  in  conjunction  with 
elastic-viscoelastic  correspondence  principle  for  the 
analysis  of  the  material  damping  for  fiber  reinforced 
composites.   Damping  values  obtained  through  these  two 
approaches  are  close  to  each  other  under  certain  situa- 
tions.  A  discontinuous  fiber  composites  model  is 
developed  to  determine  the  longitudinal  modulus  of  dis- 
continuous fiber  reinforced  composites.   Aligned  or  off- 
axis  unidirectional  fiber  composites,  in-plane  randomly 
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oriented  fiber  composites,  and  several  kinds  of  laminated 
fiber  composites  are  considered  in  damping  and  stiffness 
analysis.   Experimental  results  on  damping  and  stiffness 
by  the  impulse  hammer  technique  agree  with  the  analytical 
results  for  unidirectional  fiber  composites  and  for  cer- 
tain cross  ply  fiber  composites.   In  the  energy  approach, 
a  three-dimensional  finite  element  method  based  on  three- 
dimensional  elasticity  is  applied  to  determine  the  strain 
field  of  an  elastic  body.   Dissipated  energy  can  be  deter- 
mined through  this  strain  field  and  loss  moduli.   Damping 
through  energy  approach  depends  on  the  boundary  condi- 
tions, the  loading  conditions,  and  the  geometry  (espe- 
cially the  dimension  in  thickness  direction)  of  the  body. 

Sequential  Simplex  method,  laminated  plate  theory, 
and  an  elastic-viscoelastic  correspondence  principle  are 
used  to  optimize  both  material  damping  and  specific  stiff- 
ness of  composites.   Minimum  flexure  deformation  design  of 
orthotropic  laminated  fiber  composite  plates  is  obtainable 
through  this  optimization  procedure. 
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CHAPTER  1 
INTRODUCTION 


Due  to  the  high  strength-to-weight  and  stiffness-to- 
weight  ratios,  composite  materials  are  ideal  for  weight- 
sensitive  structures  such  as  aircraft,  spacecraft,  and 
automotive  vehicles.   In  recent  years,  with  the  advent  of 
jet  propulsion,  particularly  with  the  current  increased 
interest  in  short  take-off  and  landing  aircraft,  it  has 
become  necessary  to  pay  increasing  attention  to  the 
higher  frequency  motions  of  such  structures.   These 
motions  depend  strongly  on  the  structure's  damping  or 
capability  for  dissipation  of  vibratory  energy  [1].   In 
addition,  a  new  type  of  excitation  has  become  more  preva- 
lent, random  excitation  either  of  mechanical  or  acoustical 
origin  [2],   For  example,  jet  engine  exhaust  generally 
contains  a  noise  spectrum  wide  enough  to  excite  most  of 
the  natural  frequencies  encountered  in  aircraft  structures 
[3].   The  natural  resonance  phenomena,  so  produced,  can  be 
very  destructive.   Since  near  resonant  conditions  can  no 
longer  be  avoided  in  many  types  of  structures,  the  maximi- 
zation of  damping  within  a  structural  system  provides  a 
most  useful  concept  in  controlling  resonance  [4]. 
Unfortunately,  as  will  be  shown  in  later  chapters,  high 
damping  is  mostly  coupled  with  low  stiffness,  and  high 
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stiffness  is  mostly  coupled  with  low  damping.   Therefore, 
the  optimization  of  damping  and  of  the  stiffness-to-weight 
ratio  is  a  practical  idea  in  designing  a  proper  composite 
material  to  be  used  in  aircraft  and  space  vehicles. 

1.1   Literature  Survey 

Most  results  from  a  series  of  researches  on  damping 
beginning  in  the  1920's  [5]  indicated  that  damping  is  a 
material  property. 

Kimball  and  Lovell  [5]  experimentally  showed  that 
for  stress  cycles  of  frequency  of  from  two  to  three  a 
minute,  up  to  fifty  a  second  the  frictional  loss  (a  kind  of 
energy  loss)  is  independent  of  the  frequency  but  is  depen- 
dent on  the  amplitude  of  strain  of  the  cycles  for  eighteen 
different  solids,  including  several  metals,  glass,  cellu- 
loid, rubber  and  maple  wood,  when  strain  was  below  the" 
elastic  limit.   Crandall  [6]  pointed  out  that  the  values  for 
material  damping  which  will  be  introduced  in  Section  2.3.3 
encountered  in  practice  ranged  from  about  0.00001  to  0.2; 
however,  Lazan  [4]  pointed  out  the  material  damping  ranged 
from  0.001  to  0.1.   And  the  material  damping  depends  on 
both  the  amplitude  and  frequency  of  the  oscillation.   If, 
however,  the  system  is  completely  linear,  then  damping  is 
independent  of  amplitude  [6]. 

Recently,  the  composite  materials  got  more  attention 
in  industrial  application.   Lazan  [4]  gave  a  detailed 
review  on  material  damping  of  materials  and  material 


composites.   Kume,  Hashimoto  and  Maeda  [7]  used  the  damping- 
stress  function,  derived  by  Lazan  [8],  to  calculate  the 
material  damping  of  cantilever  beams.   They  found  that  low 
order  modes  of  a  cantilever  beam  with  equal  maximum  stress 
amplitude  gave  almost  the  same  material  damping,  theoreti- 
cally, and  experimental  results  have  the  same  order  of 
magnitude  as  the  theoretical  results  when  the  maximum  stress 
amplitude  is  less  than  a  certain  value.   Schultz  and  Tsai 
[9]  indicated  that  unidirectional  glass  fiber  reinforced 
composites  exhibit  anisotropic,  linear  viscoelastic 
behavior  when  those  undergoing  small  oscillation  and  that 
damping  increases  in  magnitude  with  change  in  fiber  orien- 
tation angle  with  respect  to  loading  direction  in  the 
order  0°,  22.5°,  90°,  45°.   Ni,  Lin,  and  Adams  [10,  11] 
used  the  laminated  plate  theory  and  two-dimensional  energy 
approach  to  predict  the  flexural  damping  of  laminated 
composites.   In  their  work,  damping  coefficient  was  deter- 
mined by  free-free  flexural  modes  of  vibration  [12].   Siu 
and  Bert  [13]  discussed  the  vibration  of  composite  plates 
having  material  damping.   Suarez,  Gibson,  and  Deobald  [14] 
observed  the  dependence  of  damping  on  frequency  of  fiber 
reinforced  epoxy  or  polyester.   Gibson  and  Plunkett  [15] 
found  that  for  small  strain,  damping  and  stiffness  are 
independent  of  amplitude  of  strains,  but,  once  the  thres- 
hold strain  is  exceeded  (i.e.,  failure  starts),  the 
resulting  increase  in  damping  is  much  more  significant 
than  the  corresponding  reduction  in  stiffness.   Similar 


results  were  also  observed  by  Tauchert  and  Hsu  [16].   Bert 
and  Clary  [17,  18]  gave  a  complete  review  on  measurement 
and  analysis  of  damping  and  dynamic  stiffness  for  compo- 
sites. The  first  paper  to  optimize  the  damping  of  the 
structure  was  perhaps  that  of  Plunkett  and  Lee  [19].   The 
damping  of  a  beam  is  improved  by  introducing  thin  con- 
strained visco-elastic  layers  on  the  top  and  bottom  of  the 
structure.   These  viscoelastic  layers  are  then  stiffened 
by  properly  designed  constraining  layers. 

Cox  [20]  discussed  the  stress  distribution  in  fibrous 
materials.   Cox's  shear  lag  stress  analysis  was  later  on 
used  to  analyze  the  stress  distribution  of  short-fiber 
composites,  as  in  studies  [21,  22].   Photoelasticity  [23, 
24]  and  finite  element  methods  [25,  26]  were  used  to 
investigate  the  stress  concentration  in  the  matrix  around 
fiber  tips  of  short-fiber  composites.    Strength  of  short- 
fiber  composites  was  analyzed  in  several  studies  [27,  28]. 
Analysis  of  complex  moduli  for  such  kind  of  material  was 
presented  in  studies  [29,  30,  31].   High  damping  of  short- 
fiber  composites  was  analytically  and  experimentally 
observed  in  references  [31,  32].   Material  damping  of 
randomly  oriented  and  unidirectional  laminar  short-fiber 
composites  has  been  discussed  by  Sun,  Wu,  Chaturvedi,  and 
Gibson  [33,  34]. 

1.2   Scope  of  This  Study 
The  objectives  of  this  study  are  to  analyze  the  mate- 
rial damping  and  to  optimize  the  specific  stiffness  (the 


ratio  of  the  stiffness  to  the  density)  and  material 
damping  of  continous  and/or  discontinuous  fibers  rein- 
forced laminated  composite  structure  elements.   The  work 
involved  in  this  research  is  briefly  introduced  as 
follows: 

A.  To  develop  a  short-fiber  composite  model  to  determine 
the  moduli  of  short-fiber  composites. 

B.  To  analyze  the  stiffness  and  material  damping  of 
unidirectional  laminar  fiber  composites,  randomly- 
oriented  fiber  composites,  and  certain  kinds  of 
laminated  fiber  composites  through  classical  lami- 
nated plate  theory  approach. 

C.  To  analyze  the  material  damping  of  laminated  fiber 
composites  by  an  energy  approach,  where  a  three- 
dimensional  displacement  finite  element  method  is 
used. 

D.  To  optimize  material  damping  and  the  specific  stiff- 
ness of  laminated  composite  plates. 

1.3   Material  Constants  and  Ranges  of  Design  Parameters 
In  this  study,  four  different  kinds  of  widely-used 
fiber  composites  (i.e.,  glass  epoxy,  Kevlar  epoxy, 
graphite  epoxy,  and  boron  epoxy)  are  involved;  and  much 
interest  is  concentrated  on  graphite  epoxy  and  Kevlar 
epoxy,  because,  generally  speaking,  the  former  has  higher 
stiffness,  while  the  latter  has  higher  damping. 


In  order  to  compare  with  the  experimental  results  of 
Suarez,  etc.  [35,  36],  the  material  constants  used  in  this 
study  are  the  same  as  those  of  the  experimental  specimens. 
Some  material  constants,  which  are  not  given  in  those 
experimental  data,  are  obtained  from  reference  [37]. 
Unless  specially  specified,  the  material  constants  used  in 
this  study  are  given  in  Table  1.1. 

The  length  of  fiber  is  one  of  the  characteristic 
parameters  of  short-fiber  composites.   Due  to  the 
existence  of  a  critical  fiber  length,  sc  (i.e.  the  minimum 
fiber  length  in  which  the  ultimate  strength  Ofu  can  be 
achieved  [38]),  there  is  a  minimum  value  for  fiber  length, 


s      °c  -i      1 

s  ->   c  _   fu  1.1 


=  > 


2Ty 


where  x  is  the  matrix  yield  stress  in  shear,  and  d  is  the 
fiber  diameter.   The  lowest  cjfu  of  those  four  fiber  com- 
posites is  2750  MPa,  as  given  in  reference  [37],  for 
graphite  T-300.   And  the  matrix  yield  shear  stress  is, 
according  to  the  manufacturer's  test  results,  97  MPa  of 
AS4/3501  graphite-epoxy  tape.   Therefore,  the  fiber  aspect 
ratio,  the  ratio  of  fiber  length  to  its  diameter,  should 
be 

s/d  >  14.2  1.2 


In  this  study,  the  fiber  aspect  ratio  is  chosen  to  be 
between  25  and  10000;  while  the  fiber  volume  fraction  V* 
is  chosen  to  be  0.65  for  the  most  cases  or  0.5  for 
randomly  oriented  fiber  composites. 


Table  1.1:   Material  Properties  Data  of  the  Matrix  and 
the  Fibers 


Constants 

Matrix 

Fibe 

rs 

Epoxy 

Glass 

Kevlar 

Graphite 

Boron 

E1    (Gpa) 

3.94 

72.4 

99.8 

175.8 

381.9 

E2  (Gpa) 

3.94 

13.8 

6.9 

13.8 

35.0 

G12 (Gpa) 

1.465 

27.6 

13.8 

27.6 

70.0 

v12 

0.345 

0.22 

0.376 

0.16 

0.21 

k  (Gpa) 

4.236 

— 

— 

— 

— 

Hi 

0.015 

0.0015 

0.C11 

0.0015 

0.0019 

P(kq/m3) 

1220. 

2539. 

1479. 

1760. 

2481. 

Si 

— 

2 

2 

2 

2 

h 

— 

1 

1 

1 

1 

where  £]_  and  £5  are  constants  used  in  Halpin-Tsai 
Equation,  n^  is  the  longitudinal  damping,  k  is  the  bulk 
modulus,  and  p  is  the  density. 


CHAPTER  2 
DAMPING 


2.1   Definition  of  Damping 
The  process  by  which  vibration  steadily  diminishes  in 
amplitude  is  called  damping.   In  many  ways,  the  assumption 
that  systems  possess  no  damping  is  a  mathematical  conven- 
ience, rather  than  a  reflection  of  physical  evidence.   In 
fact,  if  a  system  is  set  in  motion  and  allowed  to  vibrate 
freely,  the  vibration  will  eventually  die  out;  the  rate  of 
decay  depends  on  the  amount  of  damping.  This  reduction  in 
vibrating  amplitude  occurs  because  the  energy  of  the 
vibrating  system  is  dissipated  as  friction  or  heat  or  is 
transmitted  as  sound  [39]  .  And  this  is  why  damping  is  also 
interpreted  as  any  phenomenon  within  the  body  of  the 
material  where  energy  is  dissipated  [40]. 

The  concept  of  an  undamped  system  serves  not  only  a 
useful  purpose  in  analysis,  but  can  also  be  justified  in 
certain  circumstances.   For  example,  if  the  damping  is 
small  and  one  is  interested  in  the  free  vibration  of  a 
system  over  a  short  interval  time,  there  may  not  be 
sufficient  time  for  the  effect  of  damping  to  become 
noticeable.   Similarly,  for  small  damping  one  may  not  be 
able  to  notice  the  effect  of  damping  in  the  case  of  a 
system  with  harmonic  excitation,  provided  the  driving 
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frequency  is  not  in  the  neighborhood  of  any  of  the  natural 
frequencies  of  the  system  [41].   On  the  other  hand,  damp- 
ing of  a  given  system  should  be  considered  if  this  system 
is  subjected  to  vibration  near  its  resonant  frequencies 
because  damping  has  a  large  influence  on  the  amplitude  in 
the  frequency  region  near  resonance  [42]. 

2.2   Damping  Mechanism 
There  are  many  mathematical  models  representing  damp- 
ing. The  mechanism  of  damping  can  take  any  of  these 
forms  and  often  more  than  one  form  may  be  present  at  a 
time.   Therefore,  in  order  to  analyze  or  predict  the 
damping  of  a  given  system,  one  ideally  should  take  into 
account  all  possible  damping  mechanisms;  fortunately,  in 
most  practical  cases,  one  or  two  mechanisms  predominate  so 
that  one  may  neglect  the  effect  of  all  others.   Three  • 
widely  used  mathematical  models  for  damping  are  introduced 
below: 

2.2.1   Viscous  Damping 

The  viscous  damping  force  is  defined  as 


Fd  =  -CX  2.1 


where  the  constant  C  (of  dimension  force  per  unit  velocity) 
is  called  the  coefficient  of  viscous  damping.   This  type  of 
damping  occurs  in  lubricated  sliding  surfaces,  dash-pots, 
hydrolic  shock-absorbers  [43].   The  minus  sign  indicates 
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that  this  damping  force  is  always  opposite  to  the  direc- 
tion of  the  motion.   The  work  done  by  damping  force, 
namely,  dissipated  energy  during  one  cycle  of 
harmonic  motion,  x  =  A0sinwt,  will  be 


2tt 


f        rw    dx2        2 
<Vcyc  =]Fd  dx  =Jo  C  (^)  dt  =  ttCAo0)  2.2 


where  w  is  the  circular  frequency  in  radians  per  unit  time. 
Apparently,  dissipated  energy  due  to  viscous  damping  in  a 
cyclic  motion  is  proportional  to  the  frequency  and  to  the 
square  of  the  amplitude  of  motion. 

2.2.2   Dry  or  Coulomb  Damping 

This  type  of  damping  occurs  in  the  sliding  of  dry 
surfaces.  The  damping  force  during  motion  is  constant' and 
is  given,  according  to  Coulomb's  law,  by 

Fd  =  -fN  2.3 

where  N  is  the  normal  component  of  the  force  upon  the  sur- 
face of  contact,  and  f  is  the  coefficient  of  dry  friction. 
Damping  induced  by  the  joints  is  mainly  because  of  dry 
damping.   And  it  is  known  that  damping  of  built-up 
structures  (i.e.  structure  made  by  joining  together  skins, 
strings,  frames,  etc.)  could  further  be  caused  by  the 
effect  of  the  joint  [1]. 
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2.2.3   Material  Damping 

This  kind  of  damping  is  also  referred  to  as  internal, 
hysteresis  or  structual  damping.   It  is  caused  by  the 
internal  friction,  the  viscoelastic  behavior  of  the 
material,  and  the  interfacial  slip  in  the  material  itself. 

It  is  well  known  that  an  elastic  body  which  is 
repeatedly  stressed  becomes  hot.   If  an  elastic  body  is 
subjected  to  forced  oscillation,  a  positive  work  of  the 
exciting  force  must  be  spent  to  keep  the  amplitude  of  the 
oscillations  constant  in  time.   The  reason  for  the  heating 
of  the  body  and  the  expenditure  of  external  work  is  the 
internal  friction  of  the  material.   Although  this 
explanation  can  be  easily  accepted  in  a  qualitative  way, 
it  is  more  difficult  to  translate  the  problem  into 
mathematical  terms. 

Three  principal  hypotheses  have  been  proposed  to  . 
explain  the  phenomenon  of  internal  friction,  i.e.,  the  vis- 
cous theory,  the  hereditary  theory,  and  the  hysteresis  loop. 

Viscous  theory.   The  viscous  theory  assumes  that  in 
solid  bodies,  there  exist  some  viscous  actions  which  can 
be  compared  to  the  viscosity  of  fluids.   These  viscosity 
effects  are  assumed  to  be  proportional  to  the  first 
time  derivative  of  strain.   The  coefficient  of  the 
proportionality  (constant  for  each  material  at  constant 
temperature)  is  called  the  coefficient  of  viscosity. 
Based  on  this  assumption,  many  mathematical  models  (such 
as  Maxwell  model,  Kelvin-Voigt  model,  and  three-parameter 
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models    [44,    45])    have  been    introduced   to   represent  dif- 
ferent materials.      For  the   case  of  Kelvin-Voigt  model 
under  normal   deformation,    the  relationship  between  stress 
a   and  strain   e   is   expressed  by 

a  =  Ee   +   K  f|  2-4 

where  E  is  Young's  modulus  and  £  is  the  coefficient  of 
viscosity  of  the  material. 

Hereditary  theory.   The  hereditary  theory  attributes 
the  dissipation  of  energy  due  to  material  damping  to  the 
elastic  delay  by  which  the  deformation  lags  behind  the 
applied  force  [43].   According  to  this  theory,  the  defor- 
mation at  a  given  instant,  instead  of  depending  only  on  the 
actual  applied  stress  at  that  time  as  it  would  if  the 
materials  followed  Hooke's  law,  depends  on  all  the 
stresses  previously  applied  to  the  elastic  body.   The 
stress-strain  relationship  is  given  by 


ft 
a  =  Ee  +  I    <Mt,  t)  e(x)  dT  2.5 


where  t  is  the  actural  time  and  x  is  an  instant  of  time 
between   t  =  -°°  and  t  =   t.      The  function  *(t,    x)    is 
called   the   hereditary  kernel   or  memory  function. 

Hysteresis    loop.      For   a   material   under   a   cyclic 
loading,    the   stress-strain   curve   is   a   closed   curve  which 
is   called   the   hysteresis   loop.      The  physical   meaning  of 
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this  hysteresis  loop  is  given  in  Section  2.3.   The  area 
within  the  hysteresis  loop  is  proportional  to  the 
dissipated  energy.   This  area,  being  a  material  property, 
may  or  may  not  depend  on  the  frequency.   A  mathematical 
model  can  be  used  to  explain  the  energy  dissipation,  when 
this  energy  dissipation  is  independent  of  frequency  in  a 
material  loaded  by  a  cyclic  force.   In  this  mathematical 
model,  the  damping  force  is  assumed  to  be  proportional  to 
velocity  and  inversely  proportional  to  frequency,  i.e. 


F    h  x  2*6 

a   w 


If  an  external  force  Fe  is  applied  just  enough  to 
balance  the  damping  force  and  to  maintain  a  simple 
harmonic  motion,  x  =  A0sinu;t,  then 


Fe  =  Kx  +  |  i  2-7 


where  Kx  represents  the  elastic  force  of  the  system;  for 
example,  a  single  spring,  u  is  the  circular  frequency, 
A  is  the  amplitude,  and  h  is  the  hysteretic  damping 
constant.   The  relation-  ship  between  Fe  and  x  is  given  in 
Equation  2.7,  and  the  plot  of  external  force  Fe  as  a 
function  of  displacement  is  a  skewed  ellipse,  as  in 
Figure  2.1. 


(2,2  +  (FP  -  KX)2  .  x  2.8 

A0        hAQ 
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Figure  2.1:   Sketch  of  Hysteresis  Loop 

The  work  done  by  damping  force  (dissipated  energy)  in  one 
cycle  (UD)cyc  is 


(Vcyc 


Fd  dx  = 


2tt 


h     '     0  2 

-UP  dt  =  tt  h  Aq 


2.9 


Thus,  the  energy  dissipated  in  one  cycle  is  proportional 
only  to  the  square  of  the  amplitude.  This  expression 
agrees  with  the  results  of  experiments  of  Kimball  and 
Lovell  [5]  which  indicate  that  for  a  large  variety  of 
materials  such  as  metals,  glass,  rubber  and  maple  wood, 
subject  to  cyclic  stress  such  that  the  strains  remain 
below  the  elastic  limit,  the  internal  friction  is  entirely 
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dependent  on  the  rate  of  strain.   Equation  2.9  also  agrees 
with  Lazan's  notes  [46]  on  dissipated  energy.   For  example, 
at  low  amplitudes  of  stress,  the  dissipated  energy  is 
proportional  to  the  square  of  the  stress  amplitude,  and  the 
hysterestic  loop  is  elliptical  in  form. 

Unlike  homogeneous  materials,  fiber  reinforced  mate- 
rials have  interfaces  between  matrix  and  fiber.   When  a 
fiber  reinforced  composite  is  subject  to  a  tensile  strain 
cycle,  the  high  shear  stress  may  cause  the  fiber  matrix 
interface  to  fail  so  that  energy  is  dissipated  by 
friction  as  the  matrix  slides  over  the  fibers  [47].   Damp- 
ing is  then  increased  due  to  interfacial  slips  between 
matrix  and  fiber.   In  this  study,  perfect  bonding  between 
matrix  and  fiber  is  assumed;  conseqeuntly,  interfacial  slip 
is  not  considered. 

2.3   Types  of  Damping  Representations 
Many  different  disciplines  have  been  concerned  with 
damping  measurements,  and  this  has  further  complicated 
nomenclature.   Confusion  has  been  caused  not  only  by  the 
large  variety  of  damping  units  used,  but  also  by  the  lack 
of  unique  definition  for  many  well-accepted  units.   It  is, 
therefore,  desirable  to  review  the  various  damping  units 
currently  used  and  to  indicate  relationships  between  them. 
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2.3.1   Damping  Ratio  (5) 

Figure  2.2  shows  a  single  degree  of  freedom  system 
with  viscous  damping,  excited  by  force  F(t). 


Oj  c 


Pt 


F(t) 


Figure  2.2:   Sketch  of  Viscous  Damping  Model 
Its  differential  equation  of  motion  is  found  to  be 


Mx  +  Cx  +  Kx  =  F(t) 


2.10 


If  F(t)=0,  one  has  the  homogeneous  differential  equation 
whose  solution  corresponds  physically  to  that  of  free- 
damped  vibration.   The  general  solution  to  this  homogeneous 
equation  is 


-(c/2m)t  ,   bt      -bt 
x  -  e  (cxe   +  c2e   ) 


2.11 


where 


2    ,  ,1/2 
b  =  [ (c/2m)  -  k/m] 


2.12 
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cx  and  c2  are  constants  to  be  determined  by  initial 
conditions.      In  order   to  have  oscillation,    one  will   expect 


to  have 


2    <    Win  2.13 


(c/2m)      <   k/m 


Apparently,  there  exists  a  critical  value  cc  for  c, 
when  (c/2m)  equals  k/m.  Damping  ratio  [38],  £,  is 
defined  as 


where 


c  2.14 


cc  =  (4mk)  -  2.15 


2.3.2   Logarithmic  Decrement  (6) 

A  convenient  way  to  determine  the  amount  of  damping 
present  in  a  system  is  to  measure  the  rate  of  decay  of 
free  oscillations.  Logarithmic  decrement  [38]  is  defined 
as  the  natural  logarithm  of  the  ratio  of  any  two 
successive  amplitudes,  as  in  Figure  2.3,  of  a  free 
vibration. 
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Figure  2.3:   Sketch  of  Logarithmic  Decrement 

2.3.3   Loss  Tangent  (tan  $) 

It  is  well  known  that  polymer  behaves  as  viscoelastic 
material,  i.e.,  combining  two  material  properties,  one  of 
which  is  perfectly  elastic,  while  the  second  is  viscous 
fluid  [43].   Let  such  a  viscoelastic  material  be  subject 
to  a  sinusoidal  stress  experiment  at  frequency  u  such 
that  the  period  2ir/u  of  oscillation  is  sufficiently  large 
as  compared  to  the  transit  time  of  elastic  waves  through 
the  specimen  that  stress  and  strain  can  be  considered 
uniform  throughout  the  test  section.   Under  these  condi- 
tions, the  response  to  a  steady-state  sinusoidal  stress  a 
is  a  steady-state  sinusoidal  strain  e  at  the  same 
frequency  [44],  out  of  phase  by  the  angle  ♦,  e.g. 


a  =   aQ   sinut 


2.17 
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e  =  eQ  sin(a)t  -  <f>) 


2.18 


Both  the  response  amplitude  and  the  phase-shift  (or  phase 
angle)  *  are  frequency-dependent,  but  in  the  linear  range 
eQ   is  proportional  to  a  .   The  phase  relationships  are 
conveniently  shown  in  the  rotating-vector  representation  of 
simple  harmonic  motion,  as  in  Figure  2.4. 


E"e 


C.  £ 


B 


Figure  2.4 


Rotating-Vector  Representation  of 
Harmonic  Motion 


The  rotating  vector  OB  of  magnitude  E'eo  lags  behind 
the  stress  OA  by  $  radians.   Stress  OA  may  be  resolved 
into  two  components,  E'e   in  phase  with  strain  and  E"e  , 
it/2  radians  out  of  phase  with  strain,  as  in  Figure  2.4. 
Here  E'  is  the  storage  modulus  and  E"  is  the  loss  modulus. 
The  loss  tangent  tan  <J>  is  defined  as 


tan  <J>  =  E"/E 


2.19 


Some  authors  call  tan  $  the  loss  coefficient.   The  ratio 
E"/E'  is  a  measure  of  the  ratio  energy  loss  to  energy 
stored,  as  will  be  shown  in  Appendix  A.   For  viscoelastic 
material,  the  moduli  are  often  expressed  in  terms  of  a 


20 

complex  number,  called  the  complex  modulus.   In  this 
study,  the  superscript  *  is  used  to  indicate  the  complex 
modulus,  for  example: 

E*  =  E'  +  iE"  2.20 

where  i  is  J-l. 

2.3.4   Specific  Damping  Capacity  for  Cyclic  Loading  (ipc) 

The  physical  meaning  of  the  hysteresis  loop  of 
Section  2.2 .3  is  considered  here.   Since  materials  do  not 
behave  in  a  perfectly  elastic  manner  even  at  very  low 
stress  [46],  inelasticity  is  always  present  under  all 
types  of  loading,  although  in  many  cases  extremely  precise 
measurements  are  necessary  to  detect  it.   Under  a  cyclic 
loading  condition,  inelastic  behaviors  lead  to  energy  dis- 
sipation.  This  means  that  the  stress-strain  (or  load- 
deformation)  curve  is  not  a  single-valued  function  but 
forms  a  hysteresis  loop.   Energy  is  absorbed  by  the  mate- 
rial system  under  cyclic  load,  and  the  energy  absorbed  is 
proportional  to  the  area  within  the  hysteresis  loop  [46], 
as  in  Figure  2.5. 
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Figure  2.5:   Hysteresis  Loop  of  An  Inelastic  Body- 


Consequently,  another  measurement  of  damping  called 
specific  damping  capacity  [46]  can  be  obtained  by  compar- 
ing the  energy  dissipated  (or  absorbed)  (UD)cyc  of  the 
system  in  a  cycle  with  the  maximum  strain  energy  stored 


<us>max 


in  the  system  during  that  cycle. 


*c  =  (uD,cyc/(U8,«ax 


2.21 


In  Appendix  A,  viscoelastic  material  is  shown  to  have  such 
a  hysteresis  loop  under  cyclic  loading,  and  the  same 
expression  for  specific  damping  capacity  is  obtained.   The 
difference  is  that  specific  damping  of  viscoelastic 
material  is  also  a  function  of  frequency,  as  reported 
in  studies  [48,  49].   This  is  because  the  storage  and  loss 
moduli  are  functions  of  frequency. 
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For  small  damping,  the  relationships  between  those 
representations  of  damping  are  given  in  references 

[50,  51] . 
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CHAPTER  3 
DAMPING  OF  UNIDIRECTIONAL  FIBER  COMPOSITES 


3.1   Introduction 


The  objective  of  this  chapter  is  to  determine  theo- 
retically the  damping  of  unidirectional  fiber  reinforced 
polymer  matrix  composites.   The  major  damping  mechanism  of 
such  composites  is  the  viscoelastic  behavior  of  the  poly- 
mer and  fibers.   The  analysis  is  carried  out  by  first 
applying  the  concepts  of  balance  of  force  and  equal  strain 
energy  on  short-fiber  composite  model  to  determine  the 
longitudinal  modulus  of  short-fiber  composite.   Then  the 
elementary  mechanics  approach  is  used  to  find  the  modulus 
E  along  the  loading  direction  as  a  function  of  the  mecha- 
nical  properties  of  the  fiber  and  matrix  materials.   This 
is  followed  by  applying  the  viscoelastic-elastic 
correspondence  principle  [45,  52]  to  express  the  mechani- 
cal properties  of  the  composite,  fiber,  and  matrix;  then 
after  the  real  and  imaginary  parts  of  complex  modulus  are 
separated,  the  damping  of  the  composite  can  be  obtained. 
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3.2  Damping  Analysis  of  Unidirectional  Fiber  Composites 
3.2.1   Short-Fiber  Composite  Model 

The  short-fiber  composite  model  is  composed  of  a 
finite-length  fiber  and  the  polymer  matrix,  as  in 
Figure  3  .1  (a) . 
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(b) 
Figure  3.1:   Short-Fiber  Composite  Model 


Figure  3.1(c)  is  the  homogeneous  material  equivalent  to 
the  composite  of  Figure  3.1(a).   Figure  3.1(b)  is  the 
front  middle  longitudinal  section  view  of  Figure  3.1(a), 
where  d  and  s  are  the  diameter  and  the  length  of  fiber, 
respectively;  D  and  L  are  the  diameter  and  the  length  of 
the  composite  model,  respectively;  and  P  is  interpreted  as 
the  distance  between  fiber  tips  along  fiber  direction. 
The  ratio  of  P  to  s  is  defined  as  R  and  is  interpreted  as 
the  degree  of  discontinuity.   During  the  derivation  of 
Young's  modulus  along  the  fiber  direction,  the  short-fiber 
composite  model  is  treated  as  if  it  is  composed  of  two 
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materials  connectedin  series  along  the  fiber  direction. 
One  material  which  is  between  sections  H  and  H*  is  the 
mixture  of  fiber  and  matrix  having  length  s,  while  the 
other  material  is  just  the  pure  matrix  having  length  P. 

As  in  some  other  analytical  work  [23,  31,  32]  on 
short-fiber  composites,  the  results  of  Cox's  shear  lag 
stress  analysis  [20]  are  used  in  this  study.   The  expres- 
sion for  elastic  stiffness  of  the  discontinuous  fiber 
composite  is  derived  from  the  average  of  fiber  stress 
based  on  Cox's  fiber  stress  distribution  (in  which  the 
longitudinal  fiber  stress  is  a  function  of  position) . 


a.  =   CfEf    h  _  cosh[B(s/2-x)]i  3.1 

Of        £f*.f  ti     Cosh(6s/2)   ; 


where  x,  $,  and  gs/2  are  defined  in  Appendix  B,  and  ef  is 
the  strain  of  the  fiber.   In  this  study,  the  square  pack- 
ing array  of  fiber  composites  is  considered;  therefore, 
Bs/2  can  be  written,  according  to  reference  [31],  as 


G      1/2  T  ? 

£s  _  2  s  ( m    )  ->«^ 

2d    Ef  Ln  i_ 

4Vf 


The  average  fiber  stress  is 

•s/2 

af  dx  3.3 


'f 


i    r/z 
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Substitute  Equation  3.1  into  Equation  3.3 


~   _  -   r-   n  _  tanh(Bs/2)  i  3  a 

°f  "  ef  Ef  [1  ~  Bs72   ] 


For  the  composite  between  sections  H  and  H'  in  Figure 
3.1(a),  in  order  to  have  static  equilibrium,  the 
total  longitudinal  force  q  applied  to  this  composite 
must  be 


°c   Ac  =  °f   Af  +  cm  \i 


3.5 


Therefore, 


~ac   =  Ec  ec  =  af  Vf  +  am  vm  -3'6 


where  Vf  and  Vm  are  the  fiber  volume  fraction  and  matrix 
volume  fraction  within  sections  H  and  H'  separately. 
It  is  assumed  that  the  composite,  fiber,  and  matrix  (all 
between  sections  H  and  H')  have  the  same  extensional 
strain  £.   The  longitudinal  modulus  of  material  between 
sections  H  and  H'  can  be  obtained  from  Equation  3.6. 


Ec  =  Ef  v;  u  _  tanh(gs/2)]  +  Em  Vm'  3.7 
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If  Vf   and  Vm   are  expressed   in  terms  of  Vf,   Vm,    and  R,   Equa- 
tion 3.7    can  be   rewritten  as 


E      =  Ef(Vf   +   VfR)     tl    -    tanh(Bs/2)1    +   Em    (Vm  "   VfE> 

6s/2  3.8 


Alternatively,  if  the  same  assumption  is  used  as  in  conti- 
nuous fiber  composites  is  considered,  the  fiber  stress 
along  longitudinal  direction  is  assumed  to  be  uniform 
everywhere  in  fiber,  and  the  longitudinal  modulus  of 
material  between  sections  H  and  H'  can  be  obtained  by 
using  rule  of  mixtures. 

Ec  =  Ef  (Vf  +  VfR)  +  Era(Vm  -  Vf  R)  3.9 

Equations  3.8  and  3.9  show  that  for  continuous  fiber  com- 
posites, the  longitudinal  Young's  modulus  obtained  by  the 
rule  of  mixtures  is  higher  than  that  obtained  by  Cox's 
analysis.   This  is  because  in  the  rule  of  mixture 
approach,  uniform  longitudinal  fiber  stress  is  assumed, 
while  in  the  Cox's  approach,  uniform  longitudinal  fiber 
stress  exists  only  at  the  locations  far  away  from  the 
fiber  tips,  and  this  longitudinal  fiber  stress  reduces  to 
zero  at  fiber  tips.   Finite  element  stress  analyses  [25, 
26]  show  that  the  reduction  of  longitudinal  fiber  stress 
around  fiber  tips  does  exist;  and  the  magnitude  of  this 
stress  is  not  zero  but  finite.   So  it  is  hard  to  say  which 
approach  (rule  of  mixtures  or  Cox's  analysis)  is  more 
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nearly  correct.   However,  Table  3.1  shows  the  values  of 
tanh  (Bs/2) /(Bs/2)  of  graphite  epoxy  and  Kevlar  epoxy  with 
Vf  being  0.7  or  0.4.   This  table  indicates  that  the  modi- 
fication term  tanh(6s/2)  /  (Bs/2)  becomes  important  when 
fiber  volume  fraction  and  fiber  aspect  ratio  are  both 
small.   On  the  other  hand,  when  the  fiber  volume  fraction 
is  greater  than  0.4  and  the  fiber  aspect  ratio  is  greater 
than  100,  the  effect  of  tanh(6.-?/2)  (Bs/2)  could  be 
neglected. 

Table  3.1:   Influence  of  V£  and  s/d  on  tenh  (Bs/2)  (Bs/2) 


Composite 

V 

s/d 

tanh  (8: 

s/2) /(Bs/2) 

Graphite-epoxy 

0.4 

5. 

0.725 

Graphite-epoxy 

0.7 

5. 

0.369 

Kevlar-epoxy 

0.4 

5. 

0.610 

Graphite-epoxy 

0.4 

25. 

0.180 

Graphite-epoxy 

0.7 

25. 

0.074 

Kevlar-epoxy 

0.4 

25. 

0.135 

Graphite-epoxy 

0.4 

100. 

0.045 

When  the  same  external  stress  a  is  applied  to  short- 
fiber  composite  model  and  to  its  equivalent  homogeneous 
materials,  the  same  strain  energy  density  U  is  presumed 
for  those  two  materials  (Figure  3.1(a)  and  Figure  3.1(c)) 


ua  "  Uc 


3.10 
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and 


2  2 

u   =  _lo_P  +   l   £_s_  3.11 

a  '   2  E^      L      2    Ec   L 


U   = 


a2  3.12 


C      2    EL 

where  EL  is  the  longitudinal  Young's  modulus  of  the  homo- 
geneous material  equivalent  to  short-fiber  composite 
model.   After  Equations  3.11  and  3.12  are  substituted  in 
Equation  3.10,  EL  can  be  expressed  as 


3.13 


E 
c 

E 
m 

L 

E 
c 

R 

1+R 

+   E 
m 

1 
1+R 

where   Ec    can   be   obtained   by  Equation   3.8    or  Equation   3.9, 
and  R  is  the  ratio  of  P  to  s.      When  R  equals  zero,    EL 
equals  Ec.      On  the  other  hand,    if  R  is  a  very  large 
number,    Et    will   be  very  close   to  Em.      Since   the   fiber 
aspect  ratio  considered  in  this   study   is  between  25   and 
10,000,    unless   specially  mentioned,     the   modified  Cox's 
analysis    (i.e.,    Equations   3.8   and  3.13)    is  used  to 
anaylize    the    longitudinal    Young's   modulus   of    short-fiber 
composite   model. 
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It  should  be  noted  that  the  continuous  fiber  compo- 
site can  be  induced  either  by  letting  the  fiber  aspect 
ratio  be  a  very  large  number  in  the  modified  Cox's  analy- 
sis or  by  letting  R  be  zero  in  the  modified  rule  of 
mixture. 

3.2.2.   Damping  of  Aligned  Short-Fiber  Composites 

A  typical  representative  volume  element  of  off-axis, 
short-fiber  composite  is  shown  in  Figure  3.2 


Figure  3.2:   Representative  Volume  Element  of  Off-Axis 
Short-Fiber  Composites 
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For  a  continuous  aligned  composite,  the  off-axis  modulus 
Ev  along  the  loading  direction  is  given  in  reference  [38] 

A 


4       4  v       2      2 

1_  =  cos  6  +  sin  6  +  (_1 —  _  2  _LT)  sin  e  cos  0 

Ex      EL      ET       GLT       EL  _  .  . 

3  .14 


where  EL  and  ET  represent  the  moduli  along  and  transverse 
to  the  fiber  direction  respectively,  GLT  is  the  in-plane 
shear  modulus,  and  vLT  is  the  major  Poisson's  ratio. 
Equation  3.14  can  be  easily  derived  from  the  elementary 
mechanics  approach  for  off-axis  continuous  fiber 
composites.   For  off-axis  aligned  short-fiber  composites, 
one  can  derive  a  similar  expression  for  Ex  from  Equation 
3.14  by  replacing  EL,  ET,  GLT,  and  vLT  by  the  correspond- 
ing formula  for  aligned  short-fiber  composites.   The 
longitudinal  modulus  can  be  obtained  from  Equations  3.8 
and  3.13.   The  transverse  modulus  ET,  in-plane  shear 
modulus  GtT  and  the  major  Poisson's  ratio,  those  material 
constants  are  assumed  to  be  independent  of  length  of 
fiber,  can  be  obtained  by  using  the  Halpin-Tsai  Equation 
[53]  and  the  rule  of  mixtures,  i.e. 


1  +  2n  V  _,  ,_ 

E  =  E   kJ.  3'15 

T    ml-  nj^Vj 


G    =  G   IS  3'16 

LT    ml-  n2vf 
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VLT  "  vfLT  Vf  +  vm  Vm  3'17 


where 


„  -  (EfT/Em'  -  1  3.18 

1   (EfT/Em)  +  2 


n_   (GfLT/Gm'  *  1  3.19 

2    <GfLT/Gm»  +  i 


Up  to  now,  all  equations  presented  in  this  chapter 
are  derived  for  elastic  material.   When  a  viscoelastic 
material  is  considered,  the  elastic-viscoelastic 
correspondence  principle  [52]  can  be  used  to  obtain  the 
corresponding  relationships  of  viscoelastic  material.   For 
viscoelastic  material  the  basic  material  properties  are 
redefined  as 
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EfL   "  EfL  +  i  EfL 


E|T   -  E^T  +  i  E"fL 


GfLT  -  GfLT+  i  GfLT 


E*    =  E '   +  i  E"  3.20 

m      m       m 


G*    =  G'   +  i  G" 

mm       m 


v*   =  v1   +  i  v" 
mm       m 


vfLT  _  vfLT 


E*  =  E'  +  i  E"  3.21 

X      X        X 


Where  the  prime  quantities  indicate  the  storage  moduli  or 
storage  Poisson's,  the  double  prime  quantities  indicate 
the  loss  moduli  or  loss  Poisson's  ratio,  and  the  i  is 
defined  as  J-l.   In  this  research,  bulk  modulus  of  epoxy 
matrix  Km  is  assumed  to  be  real  and  independent  of  fre- 
quency [15].   For  isotropic  epoxy  matrix, 


v         -  m  3.22 

m    3(1  -  2v  ) 
m 
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While  the  viscoelastic  behavior  of  epoxy  is  considered, 


E'  +  E"  ,  „ 

_       mm  3.23 

m  "  3(1  -  2v'  -  I  2v") 
m       m 


The  complex  form  of  vm  can  be  obtained  by  Equation  3.23 


v'  +  i  v«  =   1  (1  -  — ■  ™)  3'24 

m      m    2         3  k 

m 


After  separation  of   the  real    and   imaginary  parts  of   right- 
hand   side  of   Equation   3.24,    one  will    have 


E" 

v'    +   i   v"   =   v'      +   i    [JL    (v»    -   h]  3*25 

m  mm  E '  m        2 

m 


Similarily,  the  complex  form  of  shear  modulus,  G1  +  G"   of 
viscoelastic  matrix  can  be  expressed  as  functions  of  E^, 

vm'  and  Em'  L-e-' 

E'  E'        9K       E" 

r.    +  i   R-  -    E_   +  i  m  m  m 

v    m  m      m    m  m 


m       m 


3.26 
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The  complex  form  of  Bs/2  of  Equation  3.3,  as  shown  in 
Appendix  B,  is 


B's  .  1  B"s   =  Js  ,  i    Bs  ,_m  _  _f_L}  3.27 

2        2    "   2      4   <6i    E'l; 


The  reason  that  the  imaginary  part  of  fiber  Poisson's 

ratio  is  set  as  zero  in  the  last  equation  of  Equation  3.20 

is  because  first,  most  fibers  are  known  to  be  anisotropic 

materials;  therefore,  Equation  3.22  is  not  true  for  most 

fibers.   Secondly,  the  corresponding  term  of  Em/Em  for 

most  fibers  (i.e.,  E£  /Ei)  ,  except  Kevlar  fiber,  is  much 

less  than  E"/E'  .   Due  to  the  lack  of  the  available  data 
m  m 

of  transverse  damping  and  shear  damping,  those  two  dampings 

are  assumed  to  be  equal  to  E"  /F'  . 

f  L   iL 

It  should  be  noted  that  n   and  ri   are  treated  as 

m      f 

as  material  properties,  and  they  are  defined  as 


E" 

n   =  ™  3.28 

"m   E1^ 
m 


EfL 
n  =  -4=  3.29 


By  using  Equations  3.25  and  3.28-3.29,  the  right-hand 
side  of  Equations  3.20,  3.21,  and  3.27  can  be  rewritten  as 
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E'         +    i    E"      =   E'       (1   +   i    n- 

fL  fL  fL  f 


Elm      +    i    E"      =   E'       (1   +   i    n 
fT  fT  fT  f 


GL        +    i    G"      =   G'       (1    +    i    n 

fLT  fLT         fLT  f 


E'         +    i    E"      =   E'        (1    +    i    n 
m  m  m  m 


G'         +   i   G"      =   G1       (1   +    i    n 
m  m  m  Gm 


=   E 


=   E 


=    G 


fL 

* 
fT 

* 


fLT 


=   E 


m 


=    G 


3.30 


m 


m 


+   i   v"      =   v'      +   i   n      (v' 
mm  mm 


1\  * 

i)   =  v 

2  m 


v'         =      v 
fLT  fLT 


=    v 


fLT 


E'         +      i   E"    =   E 
x  x  : 


6^+  i  JO     =|§  tl  +  i  1    (n       -  n  )]   =     |I§ 

2  2  2  2         Gm  f  2 


where  n    is  defined  as 
Gm 


9K 


IT. 


G" 

m 


'Gm 


9Km  "  Em 


nm  =  G' 


m 


3.31 


From  Equation  3.31,  it  is  observed  that  nGm  is  higher  than 
nm.  After  using  Equations  3.20,  3.27,  and  3.30,  one 
can  rewrite  Equations  3.13  and  3.15-3.17  for  viscoelastic 
material,  as  follows 
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(E'  +  i  E")  (E*  +  i  E") 
E*  =  e c m ni 3>32 

(E'  +  i  E")  ™  +  (E*  +  E")  -^ 
c      c   1+R     m    m   1+R 


1  +  2n*  V 

E*  =  (E*  +  i  E")  = — -  3.33 

T     m      m    1  -  n*  Vf 


1  +  Vf 
G    =  (G«  +  i  G")   f-1  3.34 

LT     m      m    1  -  n  V 

2  f 


V*   =  (vj   V,  +   v'  V  )  +  i  n  (v1  -  |-)V         3.35 

LT     fLT   f      mm       m   m    2   m 


where 
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E'  +  i  E"   =  (E'T  +  i  E"  )  (V,.  +  V.R) 
c      c      rL      rL   r    r 


[x  _  tanh(Bs/2)  ]    +    (E.  +  iE»  )  (v  _  v  R) 

„*  /„        m     m   m   f 
8  s/2 


tanh(6*s/2)  =  tanh(6s/2) 


+  i  |S   nGm  -  nf  1 3>37 

2  cosh   (Bs/2) 
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(E£T  +  i  E|T)/(Ei  +  i  e;>  -  1 

nl      <Ef T  +  i  EfV >  (Em  +  i  Em}  +  2  3'38 


"2  =     (GfLT  +  ^LT^m*^  +  i 


After  substitute  E£,  E^,  g£t  ,  and  v^-j   from  Equations 
3.32-3.35  for  EL,  ET,  GT,  and  vLT,  respectively,  and 

Ex+iEx  for  Ex  into  Equation  3.14,  one  obtains 


1        cos4  0     sin4  9     ,1     -  VLT.    2n   .  2n 

E'+iE"   =  TE +  9E +  (— T-  ~  2  — r)C0S  9  Sin  9 

x    x       EL  ET         GLT       EL 

3.40 


Damping  of  the  aligned  short-fiber  composite  along  x 
direction,  nx/  is  then  determined  by 


E" 
"x  =  -E^  3"41 

X 


Equations    3.40    and   3.41    show   that  material    damping   and 

stiffness  of   aligned   short-fiber  composite   are    functions 

of   material    properties  of   fiber    (i.e.,    E'   ,    E'      ,    G'      . 

fL        fT        fLT 

vfTm,    and    n   )    and  matrix    (i.e.,   E*,    v',    and   n   ),    fiber 
rjji  f  m        m  m 

aspect   ratio    (s/d) ,    fiber   volume   fraction    (Vf) ,    degree   of 
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discontinuity  (R) ,  loading  direction  (0),  and  packing 
geometry  of  fiber.   If  the  four  different  kinds  of  pre- 
packed tapes  (glass-epoxy,  Kevlar-epoxy,  graphite-epoxy, 
and  boron-epoxy)  are  used  to  make  the  fiber  composite,  the 
design  variables  utilized  to  analyze  the  stiffness  and 
material  damping  are  s/d,  0,  R,  V^,  E^  UEfL/  Ef<r'  GfLT' 
vfLT  '  vf  '  'U*'  and  E-m  Um*    Nm'  n  J>'  i'e- 


Ex  =  fl  (-f  -m'  s/d'  9'  R'  Vf}  3*42 


nx  =  f2  (Ef,  Em,  s/d,  0,  R,  Vf)  3.43 


A  similar  approach  can  be  applied  to  determine  damping 
along  y  direction,  ny,  and  stiffness  along  y  direction, 


CHAPTER  4 
DAMPING  OF  RANDOMLY  ORIENTED  SHORT-FIBER  COMPOSITES 


4.1   Introduction 
The  objective  of  this  chapter  is  to  determine 
analytically  the  material  damping  of  in-plane  randomly 
oriented  short-fiber  composites.   The  analysis  is 
carried  out  by  using  the  extension  of  the  short-fiber  com- 
posite model  and  part  of  the  results  obtained  in 
Chapter  3.   An  averaging  procedure  is  first  applied  to 
the  six  off-axis  reduced  stiffnesses  Q^    (i,  j  =  1,  2,  6) 
with  respect  to  the  angle  0  between  the  fiber  orienta- 
tion and  the  applied  load.   The  results  of  integration 
show  that  in-plane  randomly  oriented  fiber  composites 
behave  like  a  planar  isotropic  material.   By  using  the 
properties  of  isotropic  materials,  Young's  and  shear 
moduli  can  be  obtained  as  functions  of  the  reduced 
stiffnesses  Q±j  (i,  j  =  1/  2,  6).   After  the  application 
of  the  elastic-viscoelastic  correspondence  principle  and 
separation  of  the  real  and  imaginary  parts  of  the  complex 
Young's  and  shear  moduli,  material  damping  is  obtained. 
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4.2   Damping  Analysis  of  In-Plane  Randomly  Oriented 
Short-Fiber  Composites 

For  in-plane  randomly  oriented  short-fiber  compo- 
sites, no  difference  caused  by  different  direction  paral- 
lel to  the  planes  on  which  fibers  are  laid.   The  averaging 
procedure  is  one  of  the  approaches  which  will  lead  to  the 
isotrop.   Therefore,  an  averaging  by  integrating  the  six 
moduli  of  off-axis  short-fiber  composites  with  respect  to 
0  from  0=o  to  0=tt  should  be  used.   However,  from  Equation 
3.14  for  Ex  and  similar  Equations  for  Ey,  Gxy,  vx  ,  mx  and 
my  [38],  one  finds  that  it  is  not  convenient  to  integrate 
and  obtain  the  average  Ex  in  closed  form  in  terms  of 
EL,  ET,  GLT,  and  vLT.   Instead  of  integrating  the  six 
engineering  moduli,  one  can  integrate  the  six  components 
of  the  off-axis  reduced  stiffness  of  the  plane  stress 
case  Qji  (i,  j  =  1,  2,  6)  and  obtain  the  average 
Qij/  i.e. 


Qij  =  I  Qij  d0      i/  J  "  1/  2>    6         4-1 

-1  o 


The  expression  for  Q^^  (i,  j  =  1,  2,  6)  as  a  function  of  0 
can  be  found  in  reference  [38].   After  integrating  with 
respect  to  0  from  0=o  to  0=tt  and  then  dividing  each  of  the 
result  by  n,  one  obtains 
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Qll  -  Q22  -  I  <Qll  +  Q22)  +  I  (2Q66  +  Ql2>        4'2 
Q66  "  §  'Qll  +  Q22>  -  J  Ql2  +  \   Q66  4'3 

Ql2  =  \    (Qn  +  Q22>  +  !  Ql2  "  I  Q66  4'4 

Ql6  =  Q66  =  °  4-5 

It  is  easy  to  show  from  Equations  4.2-4.4  that  the  follow- 
ing relation  exists 

Ql2  +  2Q66  "  Qll  4'6 

Therefore,  after  integration,  there  are  only  two  inde- 
pendent material  constants,  namely,  Qr  and  Gr. 

Qr  =  5u  4-7 

Gr  =  Q66  4.8 

This  implies,  as  expected,  that  in-plane  randomly  oriented 
short-fiber  composites  behave  as  planar  isotropic  mate- 
rials with  two  independent  material  constants.   The 
subscript  r  represents  randomly  oriented  short-fiber 
composites. 
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For  isotropic  materials,  the  following  relations 
exist 


Qr  =  -^V  4.9 

1-v 


2L  =  S£ 4.10 


'r 


2(1  +  v) 


where   v   is  defined  by  Qi2/Qn*      Elimination  of  v  from 
Equations   4.9    and   4.10    yields    the   expression   for 
Er,    the   Young's   modulus   of   a   randomly  oriented   short- 
fiber   composite   as    a   function  of   Gr   and   Qr. 


Er   =   4    Gr    (1   -  — ^ — )  4.11 

5r 


Substitution  of  Equations  4.2,  4.3,  4.7,  and  4.8  in  Equa- 
tion 4.11  yields  Er  as  a  function  of  the  four  reduced 
stiffness  Q-q,  Q22»  Q].2'  and  Q66' 


Er  =  [§  (QX1  +  Q22)  -  Q12  +  2Q66] 


[1  - 


Qll  +  Q22  -  2Q12  +  4Q66 
3(Qn  +  Q22)  +  ^  (Q12  +  2Q66) 


4.12 


Similarly, 
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gr  =  £  (Qn  +  Q22)  -  ?Qi2  +  bee  4-13 


Since  Qii/  ^22'  ^12'  an^  ^66  are  directly  related  to  the 
four  basic  engineering  constants  EL,  ET,  GLT,  and  vLT 
(defined  in  Equations  3.13,  3.15,  3.16,  and  3.17),  accord- 
ing to  Equations  4.12  and  4.13,  Er  and  Gr  can  be  expressed 
as  functions  of  EL,  ET,  GLT,  and  vLT. 

Next,  as  in  Section  3.2.2  for  aligned  short-fiber  com- 
posites, according  to  the  elastic-viscoelastic  correspond- 
ence principle,  one  may  replace  Er  by  E*  =  E^.  +  i  E£,  §r 
by  G*  =  G'r+  iGjl'  ,  El   by  E*  ,  ET   by  E*   ,  G^   by   G*T  , 
vLT  by  v£T  where  E£  ,  E|  ,  G£T  ,  and  v£T  ,  are  defined  in 
Equations  3.32-3.35.   After  separation  of  the  real  and 
imaginary  parts,  the  material  damping  constants  nr  and  rigr 
of  in-plane  randomly  oriented  short-fiber  composites  can 
be  obtained. 


;■■ 


gr 
nr  =  — £  4.14 

gr 


'Gr 


Gr 

_£  4.15 


2r 


An  alternative  way  of  determining  Er  and  Gr  is  given  in 
Appendix  C. 


CHAPTER  5 
DAMPING  OF  LAMINATED  FIBER  COMPOSITES- 
LAMINATED  PLATE  THEORY  APPROACH 


5.1  Introduction 

In  this  study,  laminated  plate  theory  and  an  energy 
approach  are  used  to  analyze  the  material  damping  and 
stiffness  of  symmetrically  laminated  fiber  composites.   In 
this  chapter,  we  will  discuss  all  analytical  work  of 
laminated  plate  theory  approach,  while  in  Chapter  6  the 
energy  approach  will  be  presented. 

According  to  laminated  plate  theory,  the  constitutive 
equations  (Equation  5.1)  have  already  been  given  in  refer- 
ences [38,  50],  in  terms  of  [A],  [B],  and  [D]  (i.e.  [A]*, 
[B]*,  and  [D]*)  matrices.   Material  damping  of  laminated 
composites  can  then  be  derived  from  the  expression  of  [A], 
[B] ,  and  [D] . 

5.2  Damping  Analysis  of  Laminated  Fiber  Composites 
Through  Laminated  Plate  Theory  Approach 

For  a  laminated  fiber  composite  plate,  as  in  Figure 

5.1,  the  constitutive  equations  are  given  in  references  [38, 

50]  as  shown  in  Equation  5.1. 
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Figure  5.1:   Sketch  of  Laminated  Fiber  Composite  Plate 
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In  Equation  5.1,  e°,  e°,  and  y°   are  middle  plane  strains, 

*   y  aj 

kv,  k,,  and  kv„  are  plate  curvatures  defined  in  [38,  50], 

x.  y        Ay 

and  A^-,  B^^,  and  D^  (i,  j  =  1,  2,  6)  are  the  equivalent 
reduced  in-plane  stiffness,  coupling  stiffness  and 
reduced  flexure  stiffness,  respectively.   They  are  expressed 
in  terms  of  Q^^    and  total  thickness  of  the  plate,  h,  as 
follows : 


Dij 


h/2 
-h/2 


Aij    =      I  Qij    dz      *        *       <Qij>k    (hk~hk-l> 


Bij  z  Qij   dz     =       Z     2(Qlj)k    ChJ-hJ.!)- 

K     1  D  •  a  D 

-h/2 


h/2  n      1  -  ,o 


•  Qij   dz     =       I     3<Qij)k    (hj-hj.i) 

iv~~X 


-h/2 


5.2c 
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where  Q,^    can  be  expressed  in  vector  form  [50] 


*              * 

'   Ul 

u2 

Q22 

Ul 

-u2 

Q12 

>  = 

U4 

0 

Q66 

U5 

0 

Ql6 

0 

i". 

Q26 

0 

h 

and 


u- 


0-. 


-u- 


-u- 


u- 


-u, 


1     1 


cos  2  G 


cos  4  0 


U±    (i  =  1,  2  .  .  5)  are  defined  in  [50]  as 

ui  =  I   (3Qn  +  3Q22  +  2Qi2  +  4cW 


5.3 


2  =  «  (4Qn  -  4Q22) 


Uo  = 


U3  =  -  (Qn  +  Q22  -  2Q12  -  4Q66) 


5.4 


U4  "  o    (Qll  +  Q22  +  6Ql2  "  4Q66} 


U5  =  -  (Qn  +  Q22  -  2Q12  +  4Q66) 


49 


Qij  (i/  J 
50], 


=  1,  2,  6)  are  known  to  be,  by  references  [38, 


Q 


11    l-vLT  vTL 


"22    l-vLT  vTL 


)    =    VLT  ET 
•12  '   l-vLT  vTL 


5.5 


Q66  = 


=  G 


LT 


and 


v  = 
TL 


vLT  ET 


.6 


where  EL,  ET,  GLT,  and  vLT  are  defined  in  Equations  3.13, 
3.15-3.17.   Interesting  relations  should  be  noted  that 
Q.ji  (i,  j  =  1,  2,  6)  of  Equation  5.3  will  be  equal  to  Qij 
of  Equations  4.2-4.5,  providing  U2  and  U3  are  set  as  zero. 
Equation  5.1  can  be  rewritten  in  matrix  form  as 


N 


>  = 


M 


B 


B 


D 


<    > 
k 

v    J 


5.7 
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For  symmetric  laminates,  the  coupling  stiffness  matrix  B 
is  a  zero  matrix,  and  Equations  5.1  and  5.7  can  be 
uncoupled  as 


"Nx  ' 

]Ny 

>   = 

Nxy 

All    A12    A16 


A12    A22    A26 


L16 


A26    A66 


fcy 


:  > 


fxy 


5.8a 


or 


{N} 


[A]  {e°} 


5.8b 


and 
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x 


<!my 


M 


xyj 


Dll    D12    D16 


D12    D22    D26 


D16    D26    D66 


'kx' 

< 

ky  f 

kxy 

5.9a 


or 


{M}   =    [D]  {k} 


5.9b 


From  Equation  5.8a  and  the  definition  of  A.^,  one  can 
express  in-plane  moduli  E^  (i,  j  =  1,  2,  6)  as 


-1 


_1_ 
h 


5.10 
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where  h  is  the  total  thickness  of  laminated  composite  and 
A^  are  the  elements  of  the  inverse  of  matarix  [Ajj],  and 
A^3[  are  determined  by  the  following  system  of  equations. 


r  -1   -1   -1 
All  A12   A16 


A26 


A12   A22 


-1    -1   -1 
A16   A26   A66 


kll 


L12 


L12 


L22 


L16 


l26 


A16   A26    A66 


10   0 


0   10 


0   0   1 


5.11 


According  to  elastic-viscoelastic  correspondence 
principle,  when  the  viscoelastic  behaviors  of  fiber 
composites  are  considered  the  elastic  material  constants 
EL,  ET,  GLT,  and  vLT  (defined  in  Equations  3.13,  3.15-3.17) 
should  be  replaced  by  the  corresponding  complex  moduli, 
respectively.   Consequently,  the  complex  form  of  Equations 
5.8  and  5.9  can  be  written  as 


r  * 
N„ 


N„  >  = 


N 


xy 


Lll 


^12 


'16 


^12 


l22 


^26 


*   N 

L16 


l26 


l66 
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^xy 


> 


5.12a 


or 
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{N*}    =    [A*]  {€'} 


5.12b 


r     * 


<M 


>  = 


M 
I  xyj 


'11 
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{M*} 


[D  ]  {k} 


5.13b 


The  corresponding  complex  form  of  in-plane  complex  moduli 

* 
Ejj  can  be  expressed  as 


'ID 


(Aij) 


"I 


=  Eij  +  iEij 


5.14 


*   -1 
where  (A^j)    is  the  element  in  the  ith  row  and  jth  col 

of  the  inverse  of  matrix  [A^  ]  .   Appendix  D  shows  the 

*   -1 
detail  of  the  derivation  of  (A^) 

The  in-plane  material  damping  jflji  and  flexure 

material  damping  pH^-;  are  defined  as  follows: 


umn 


!niD 


•11 


E'  . 
ID 


i,  j  =  1,  2,  6 


5.15 
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N 

P11       i,  j  =  1,  2,  6  5.16 


Fnij     D' 


For  the  same  kind  of  fiber  composite,  the  result 
obtained  by  Equation  5.15  for  sixteen  plies  of  unidirec- 
tional laminated  composite  is  same  as  that  obtained  by 
Equation  3.41  for  laminar  composite  with  the  same 
off-axis  angle  0.   This  indicates  that  the  approach 
presented  in  this  chapter  is  correct,  although  it  may  not 
be  convenient  because  the  inverse  of  a  complex  matrix  is 
involved. 


CHAPTER  6 
DAMPING  OF  LAMINATED  FIBER  COMPOSITES— ENERGY  APPROACH 


6.1   Introduction 

The  drawback  of  damping  analysis  using  laminated 
plate  theory  approach  is  that  it  does  not  include  the 
effect  of  interlaminar  stresses  (the  stresses  at  the 
interfaces  of  laminated  composites).   It  has  been  shown  in 
study  [54]  that  even  when  an  in-plane  uniform  tension  load 
is  applied  to  the  laminated  composites,  there  do  exist 
appreciable  interlaminar  stresses  around  the  free  edges. 

In  this  chapter,  an  energy  approach  in  conjunction 
with  a  three-dimensional  finite-element  method  [55,  56]  is 
used  to  analyze  the  material  damping  under  certain  loading 
and  boundary  conditions.   It  is  believed  that  this  model 
represents  a  more  realistic  approach  by  including  the 
energy  dissipated  at  the  interfaces.   By  using  this 
approach,  we  improve  the  approach  from  a  two-dimensional, 
classical,  laminated  plate  theory  to  a  three-dimensional 
elasticity  theory. 


54 


55 


6.2      Damping  Analysis   of   Laminated  Fiber  Composites 
Through  Energy  Approach 

The  damping  of  laminated  materials   in  the  first  mode 

vibration  is  determined  as 

*    <qV  eye 

n  =  — 3jj 6.1 

2    2*    qUs 
q=l         4 

where  n  is  the  total  number  of  the  plies,  (qUD)  cyc  is  the 
energy  dissipated  in  the  qth  layer  during  a  cycle,  and  qUg 
is  the  maximum  strain  energy  stored  in  the  qth  layer. 

Detailed  expressions  of  (qU^cyc  and  qUs  are  9iven  in 
Appendix  E  in  which  the  energy  expressions  for  a  visco- 
elastic  material  given  in  reference  [57]  are  used. 

The  analytical  expression  of  the  maximum  strain 
energy  D*s  for  an  elastic  body  is 

Us  =  \  Cjk   ek  dv 


(j,  k  =  1,  2  .  .  .  6) 


where   e^   are  the  maximum   strains,    and  C^  are   the  moduli 
of    the    elastic    body. 

The  expressions   of  C^k    (j,    k  =  1,    2   ...    6)    for 
orthotropic  elastic  material    are   given  by  Jones    [58]; 
here  k   and  j    axes   are  material   principal   axes.      In  this 
study,    each   layer  of   fiber  composite   is   considered   to  be 
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transversely  isotropic  material.   Therefore,  as  presented 
in  Appendix  F.l,  the  moduli  of  mth  layer  fiber  composite 
could  be  defined  as  functions 


(Cjk*m  ~  f  (EL  '  ET  '  GLT  '  VLT  '  vTT'^m         6*3 


where  EL  ,  ET  ,  GLT  ,  and  vLT  are  given  in  Section  3.1  for 

short-fiber  composites,  and  vTT,  is  assumed  equal  to  vm. 

An  analytical  expression  for  vTTi  is  available  in  the 

literature  [38],  but  the  assumption  that  vTTi  =  vm,  is 

believed  to  be  accurate  enough  in  our  analysis.   According 

to  the  elastic-viscoelastic  correspondence  principle,  when 

the  viscoelastic  behaviors  of  the  body  are  considered,  EL, 

*    *   *    * 
ET'  GLT'  VLT'  VTT'  sh°ul<3  be  replaced  by  EL,  ET,  GLT,  vLT, 

* 
and  vTT',  respectively.   For  example, 


vtl  EL 

r        = ^ h 6.4 

12   1  "  vTL  vLT 


chanqes  into 


r*  (VTL  +  iVTL)  CL   +    i  hL)  ,   . 


or 


* 
C12  =  c12  +  i  C"l2  6*6 
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After  separating  the   real   and   imaginary  parts   and  using 
matrix  rotation    (if   the  material   principal   axes  do  not 
coincide  with  the  global   axes),    the   complex  moduli  C* 
with  respect  to  the  global  axes   (i.e.,   x,   y,   and  z  axes) 
of   the  kth   layer  of   fiber   composites   are   then  obtained  in 
the   form 


_* 


[c   ]k  =    [C'  ]k  +  i    [cn]k 


6.7 


where  C  are  storage  moduli,  and  C"  are  loss  moduli.   All 
three  matrices  are  symmetric  matrices. 

The  storage  energy  qUs  and  the  dissipated  energy 
during  a  cycle  (  Ug)  c  in  the  qth  layer  of  laminated 
composites  can  be  expressed  as 


qUs 


1 
2 


ek  dv 


J   Cjk 


6.8 


and 


(qVcyc  ■  * 


ej   Cjk  ek  dv 


6.9 


correspondingly.   Consequently,  the  material  damping  of 
a  n  layers  laminated  composite  in  the  first  mode  vibration 
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Z      \         {e}T  [C»]  {e}  dv 

n-  q=nl      q 6.10 


I      {e}T  [CM  {e}  dv 
q-1  Jv0 


where  strain  field  {e}  depends  on  the  loading  and  boundary 
conditions,  and  the  strain  field  would  be  determined  in 
this  study  through  a  three-dimensional  finite  element 
method. 

From  Equation  6.10,  one  may  thus  arrive  that  even  a 
highly  dissipative  modulus  cannot  contribute  significantly 
to  the  total  loss  factor  n#  if  it's  associated  strain  (or 
stress)  does  not  participate  considerably  in  the  total 
stored  energy. 

The  procedures  taken  in  energy  approach  are  briefly 
described  as  follows: 

Step  1; 

As  in  the  laminated  theory  approach,  the  elastic 
solution  is  first  sought  in  energy  approach.   The 
equation  of  motion  of  an  elastic  body  in  static  case  can 
be  derived  from  the  principle  of  stationary  potential 
energy  given  in  reference  [59] . 


6(U-We)  -  0  6.11 


where  U  and  We  are  the  strain  energy  and  the  work  done  by 
external  forces,  respectively. 
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Step  2; 

After  the  assumed  displacement  function  (as  a 
function  of  nodal  displacement  q)  is  substituted  into 
Equation  6.11,  the  equations  of  motion  can  be  expressed  as 
a  system  of  equations  (details  are  given  in  Appendix  F.2). 

[k]  {q}  =  {?}  6.12 

where  [k]  is  the  stiffness  matrix,  and  {f}  is  the  nodal 

forces. 

Step  3 ; 

After  the  displacement  field  is  determined  by  sub- 
stituting the  solution  of  Equation  6.12  into  the  assumed 
displacement  function,  the  strain  field  {e}  of  the  elastic 
body  could  be  obtained  through  displacement-strain 
relations. 
Step  4: 

Once  the  strain  field  of  an  elastic  body  is  known,  the 
material  damping  of  a  viscoelastic  body  under  the  same 
loading  and  boundary  conditions  can  be  determined  by 
Equation  6.10. 

6.3   Loading  and  Boundary  Conditions 

Two  sets  of  loading  and  boundary  conditions  are  dis- 
cussed in  this  study.   However,  other  loading  and  boundary 
conditions  can  also  be  accommodated  in  this  approach. 
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Case  1;    Uniform  in-plane  tension  load. 

In  this  case,  only  one  quarter  of  the  plate  is 
analyzed  by  finite  element  method.   External  stress  a 
is  applied  at  one  edge. 
Case  2:   One  side  clamped  plate  under  Mv  load. 

In  this  case,  the  plate  is  clamped  at  x=0.   The 
external  force  F  is  applied  along  x=L,  and  the  external 
force  -F  is  applied  along  x=0.983L. 

The  numerical  results  are  presented  in  Section  9.5, 


CHAPTER  7 
EXPERIMENTAL  MEASUREMENT  OF  DAMPING 


7.1   Introduction 

Free  vibration  decay  [9],  band-width  method  [91, 
resonant-dwell  method  [60],  forced-vibration  techniques 
[61],  and  impulse  techniques  [62,  63,  64]  are  the  very 
popular  experimental  techniques  used  to  measure  the 
material  damping.  All  of  these  tests  are  subject  to  the 
air  drag  [65,  66],  if  the  tests  are  not  conducted  under  a 
vacuum  condition.   In  this  study,  an  improved  impulse 
technique  approach  [64]  is  utilized  to  measure  the 
material  damping. 

The  impulse  technique  consists  of  the  application  of 
a  force  pulse  at  a  point  on  the  test  structure  and  the 
measurement  of  the  response  at  another  point.  The  input 
force  and  response  signal  are  digitally  processed  by  the 
analyzer  to  form  the  frequency  response  function  or 
transfer  function.   Damping  and  natural  frequencies 
can  then  be  extracted  from  the  output  of  the  analyzer. 
All  measures  are  performed  using  a  digital  signal  process- 
ing technique.   It  involves  filtering  and  sampling  the 
input  wave  forms.   The  sampling  process  converts  a 
voltage  (at  a  certain  point  in  time)  into  a  numerical 
representation. 
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These  numbers  are  then  processed  digitally  to 
produce  the  various  calculations  performed  by  the 
analyzer.  The  primary  output  from  the  Fourier  analyzer 
is  the  frequency  response  function,  which  is  a  measure  of 
system's  characteristics.   The  analyzer  also  calculates 
the  coherence  function,  which  ranges  from  0.0  to  1.0.   The 
larger  the  potential  measurement  noise  or  the  system  non- 
linearity,  the  lower  the  coherence  function  will  be  [67]. 
A  coherence  value  of  unity  indicates  that  the  output  is 
completely  related  to  the  input.   Thus,  the  coherence 
function  used  here  is  a  measure  of  the  "quality"  of  the 
data. 

7.2   Apparatus 


The  composite  plates  or  beams  are  fixed  by  two  alumi- 
num blocks  at  one  end.   A  non-contact  probe  (KD-2310-3U, 
Kaman  Science  Corporation,  Colorado  Springs,  Colorado), 
known  as  a  motion  transducer,  is  located  about  1.5  mm 
below  the  tip  of  each  specimen.   This  motion  transducer 
operates  under  the  principle  of  eddy  current.  An  aluminum 
foil  target  of  diameter  20  mm  is  cemented  underneath  the 
tip  of  the  specimen.  A  force  transducer  is  mounted  to  the 
head  of  the  impulse  hammer  (Model  K291A,  Piezotronics, 
Inc.,  New  York),  to  measure  the  force  input  to  the  speci- 
men. The  other  tip  of  the  hammer  is  connected  to  the 
fixed  spring  so  that  the  magnitude,  location,  and  the 
dwelling  time  of  the  impact  load  can  be  controlled.   To 
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improve  the  coherence,  an  octave  filter  (4302  dual  24db, 
Ithaco,  Inc.,  New  York),  which  is  connected  to  the  force 
transducer  at  the  hammer  tip  on  one  end  and  to  the  Fast 
Fourier  Analyzer  (FTT,  Model  5420,  Hewlett  Packard)  on  the 
other  end,  is  used  to  magnify  the  input  signal  100  times. 
The  impact  point  is  near  the  stif fer  portion  of  the  speci- 
men.  The  excitation  and  response  signals  are  fed  into  the 
FFT  analyzer,  which  displays  the  frequency  response  func- 
tion and  coherence  function.   Each  frequency  response  and 
each  coherence  frunction  are  based  on  a  statistical  analy- 
sis of  an  ensemble  of  six  tests.  A  schematic  drawing  of 
the  experimental  set-up  is  shown  in  Figure  7.1. 

The  frequency  response  function  is  a  complex  valued 
function.   A  typical  experimental  display  of  the  real  and 
imaginary  parts  of  the  frequency  response  function  for  a 
graphite-epoxy  composite,  according  to  reference  [64].,  is 
shown  in  Figure  7.2  and  an  enlarged  schematic  drawing  of 
real  part  of  frequency  response  function  is  shown  in 
Figure  7.3.   As  prescribed  in  references  [62,  63,  64],  the 
peak  of  the  imaginary  part  determines  the  resonant  fre- 
quency, and  then  from  the  corresponding  real  part  (see 
Figure  7.3)  the  material  damping  (loss  factor,  n)  can  be 
calculated  by 


2 


(fa/fb)2+  1 


n  =    ■  n-  h' £_  7.1 
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In  using  the  impulse  technique  to  measure  the  mate- 
rial damping,  following  precautions,  as  reported  in  study 
[64],  must  be  taken. 

1.  In  order  to  improve  the  accuracy  of  experimental  data 
of  damping  associated  with  a  particular  mode  of 
vibration,  the  values  of  natural  frequency  within  the 
frequency  range  of  interest  (say  0  to  1,600  Hz)  are 
first  approximately  determined.   Then  a  zooming 
technique  to  increase  the  resolution  of  the  response 
in  the  neighborhood  of  this  particular  frequency  is 
used. 

2.  It  is  necessary  to  avoid  measurement  of  response  near 
a  nodal  point  for  the  modes  to  be  tested.   Such  meas- 
urements would  consist  primarily  of  noise,  since  the 
actural  response  is  very  small  near  nodal  points. 

3.  The  amplitude  of  the  vibration  must  be  kept  below  the 
thickness  of  the  specimen  to  ensure  that  air  damping 
(air  drag)  is  negligible,  since  the  air  damping  is 
linearly  proportional  to  the  amplitude-to-thickness 
ratio  of  the  beam  specimen  [66] . 

4.  It  is  important  to  optimize  the  Analog  to  Digital 
Converter  (ADC)  range  setting  on  an  FFT  analyzer 
before  making  a  measurement,  since  an  optimized  ADC 
range  set  will  increase  resolution  in  the  digitizing 
process  [35]. 

Recently,  Suarez  and  Gibson  [35,  36],  did  some 
experiments  on  material  damping  of  short-fiber  composite 
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through  impulse  hammer  technique.      Some  of  those  experi- 
mental   results   for  unidirectional   short-fiber  composites 
presented  in  this  chapter  are   compared  with  the  analytical 
results. 
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Figure   7.3 


Enlarged  Schematic  Drawing  of  Real  Part  of 
Frequency  Response  Function 


CHAPTER  8 
OPTIMIZATION  OF  DAMPING  AND  SPECIFIC  STIFFNESS 
FOR  FIBER  COMPOSITES 


8.1   Introduction 

Plunkett  and  Lee  [19]  pointed  out  that  maximum  damp- 
ing can  be  obtained  through  properly  choosing  the  length 
and  spacing  of  the  constraining  layers.   Those  layers  are 
used  to  constrain  the  viscoelastic  layer  (or  coating)  on 
the  structure  surface.   In  the  study  presented  here,  the 
optimization  of  damping  is  based  on  the  composite  material 
itself  (i.e.,  choosing  the  proper  fiber  length,  fiber 
direction,  the  longitudinal  distance  between  fiber  tips, 
etc.).   It  is  known  that  high  damping  can  reduce  the  - 
displacement  at  and  near  the  resonant  frequencies,  while 
high  stiffness  can  further  reduce  the  displacement  at 
other  frequencies.   The  best  way  to  optimize  both  damping 
and  stiffness  is  to  keep  them  as  high  as  possible.   Unfor- 
tunately, there  exists  a  general  trend  that  higher  damping 
is  mostly  coupled  with  low  stiffness  and  vice  versa. 
Therefore,  the  idea  of  the  optimization  in  this  study  is 
to  try  to  increase  damping  without  sacrificing  the 
stiffness  too  much. 

From  the  formulations  presented  in  previous  chap- 
ters, one  can  see  the  complication  contained  in  damping 
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analysis.   It  will  be  a  simplier  approach  for  the  optimi- 
zation analysis  if  no  derivative  is  needed.   Consequently, 
the  so-called  Sequential  Simplex  Method  [68,  69]  has  been 
selected  to  analyze  the  optimization. 

8.2   Brief  Introduction  to  Sequential  Simplex  Method 
The  sequential  Simplex  Method  [68,  69]  takes  a  regu- 
lar geometric  figure  (known  as  a  simplex)  as  a  base. 
Thus,  in  two  dimensions  (i.e.,  two-design  variables),  one 
should  choose  an  equilateral  triangle;  and  in  three  dimen- 
sions (i.e.,  three-design  variables),  one  should  choose  a 
tetrahedron. 

Observations  (experiments)  are  located  so  that  the 
objective  function  is  evaluated  at  the  points  formed  by 
vertices  of  the  geometric  figure.   One  vertex  is  then 
rejected  as  being  inferior  in  value  to  the  others.   The 
general  direction  of  search  may  then  be  taken  in  a  direc- 
tion away  from  this  worst  point,  the  direction  being 
chosen  so  that  the  movement  passes  through  the  center  of 
gravity  of  the  remaining  points.   A  new  point  is  then 
selected  along  this  direction  so  as  to  preserve  the  geome- 
tric shape  of  the  figure,  and  the  function  is  evaluated 
anew  at  this  point.   The  method  proceeds  with  this  process 
of  vertex  rejection  and  regeneration  until  the  figure 
straddles  the  optimum.   When  no  subsequent  moves  would 
lead  to  further  improvement,  the  last  few  geometric 
figures  are  essentially  repeated.   There  are  three  rules 
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that  govern  the  whole  procedure  of  this  method.   These 
rules  are  explained  in  the  two-dimensional  case  as 
follows : 

1.  Take  the  point  to  be  rejected  where  the  worst  value 
of  the  objective  function  is  obtained,  and  replace  it 
by  its  reflection  in  the  opposite  side  of  the  tri- 
angle.  See  Figure  8.1,  where  point  B  is  replaced  by 
point  D. 

2.  No  return  can  be  made  to  points  which  have  just  been 
left;  see  Figure  8.2,  where  point  A  is  replaced  by 
point  D.   If  point  D  is  still  the  worst  point  of 
triangle  BCD,  then  point  B  is  replaced  by  point  E. 

3.  If  the  best  valued  vertex  remains  unchanged  for  more 
than  M  iterations,  then  the  simplex  size  is  reduced 
by,  for  example,  halving  the  distance  of  all  other 
vertices  from  that  vertex.   The  next  stage  can  then 
start.   For  example,  as  in  Figure  8.3,  point  C  is 
sequentially  repeated  in  five  triangles  ABC,  BDC, 
DEC,  EFC,  and  FGC;  and  after  the  size  of  the  fifth 
triangle  FGC  is  reduced  to  triangle  F'G'C,  the  rule  1 
is  applied  at  the  new  triangel  F'G'C. 
The  magnitude  of  M  depends  on  the  number  of 

variables.   Spendley  [68]  suggested  that 


M  =  1.65n  +  0.05  n2  8.1 


72 

where  n  is  the  total  number  of  design  variables.  The 
search  can  finally  be  stopped  when  the  simplex  is  small 
enough  to  locate  the  optimum  adequately. 

It  should  be  noted  that  Rule  1  and  Rule  2  force  the 
simplexes  to  circle  continuously  about  an  indicated  opti- 
mum, rather  than  oscillate  over  a  limited  range  such  as  a 
ridge  (see  Figure  8.2).   If  one  ensures  that  at  any  point 
violating  a  constraint,  a  sufficiently  large  negative 
value  (or  positive  value,  if  one  is  maximizing)  is  set  for 
the  response  at  such  a  point,  the  system  of  simplexes  will 
then  move  along  rather  than  cross  the  constraints. 

The  advantages  of  this  method  are  as  follows: 

1.  It  is  easy  to  apply. 

2.  It  is  useful  when  analytical  or  numerical 
derivatives  of  the  objective  function  are  not 
available. 

The  disadvantage  of  this  method  is  that  it  could  not 
guarantee  that  the  global  optimum  design  is  obtained. 
Therefore,  several  different  initial  locations  (designs) 
should  be  considered  to  get  global  optimization. 
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Xi 


Figure  8.1:   Rule  1  of  Sequential  Simplex  method 


x. 


Figure  8.2:   Rule  2  of  Sequential  Simplex  method 


X, 
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X-, 


Figure  8.3:   Rule  3  of  Sequential  Simplex  method 


8.3  Mathematical  Formulation  of  Design  Problems 
Energy  approach  and  laminated  plate  theory  approach 
are  both  utilized  in  damping  analysis.  The  main  advantage 
of  the  energy  approach  over  the  laminated  plate  theory 
approach  is  that  the  interlaminar  stresses  are  included, 
The  disadvantage  of  the  energy  approach  is  that  it  takes 
much  more  computer  time  than  the  laminated  plate  theory 
approach  does.   The  influence  of  interlaminar  stresses 
can  be  neglected,  if  the  thickness  of  each  group  of  lami- 
nated composite  is  thin  (a  detailed  discussion  is  given 
in  Section  9.5).   Then  laminated  plate  theory  approach 
for  damping  analysis  is  a  good  way  for  the  optimization 
analysis  for  thin  structures. 
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Generally  speaking,  the  design  variables  of  continu- 
ous or  discontinuous  fiber  composites  include  fiber  volume 
fraction  Vf,  fiber  moduli  Ef,  fiber  aspect  ratio  s/d,  the 
longitudinal  distance  between  fiber  tips  p,  and  fiber 
orientation  angle  0  for  certain  special  stacking  sequence. 
Among  those  six  variables,  only  fiber  aspect  ratio  s/d  and 
fiber  tip  longitudinal  distance  p  are  treated  as  the 
design  variables  in  the  optimization  analysis.   Since 
orthotropic  materials  are  utilized  in  the  optimization 
analysis,  the  fiber  direction  0  is  kept  as  certain  con- 
stants.  The  fiber  volume  fraction  is  also  set  as  a 
constant  (0.65)  in  this  study,  since  the  fiber  composites 
are  made  from  prepreg  tape.   Stacking  sequence  is  fixed 
and  most  of  this  study  is  concentrated  on  graphite-epoxy, 
Kevlar-epoxy  and  the  hybrid  of  those  two  composite 
materials,  so  that  only  two  sets  of  fiber  moduli  are  . 
considered. 

Two  different  cases  are  considered  in  optimization 
analysis: 

Case  one  is  to  optimize  the  specific  stiffness  and 
material  damping  of  an  orthotropic  square  plate  simply 
supported  on  all  four  sides  under  free  vibration. 

The  equation  of  motion  for  this  case  is 


Dll  W,xxxx  +  2<D12  +  2D66>  W,xxyy  +  D22w,yyyy  =  Pw, tt 

8.2 
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where  the  D^  are  derived  from  laminated  plate  theory,  and 
W  is  the  flexural  deformation.   Let 

w  =  Vi  Fm  <x>  Hn  Cy>  T(fc)  8'3 

By  separation  of  variables,  the  form  (8.4)  is  obtained  for 
T(t);  while  Fm  (x)  and  Hn  (y)  satisfy  Equation  8.5 


T(t)  =  e  (i  JL.)'  t     k  =  1,  2,  3  .  .  .  . 

p  8.4 


Dll  Fm  (iV)Rn  +  2<D12  +  2D66>  FmHn  +  D22  Vi^' 

=  XkFmHn  Q  _ 

O  .  D 


where  Xk  is  the  constant  to  be  determined,  and  i  is  j-1 
in  Equation  8.4. 

In  accordance  with  the  simple  supoort  boundary 
conditions,  Fm  and  Hn  functions  can  be  assumed  as 


,-.  /  \     •    miry  Q  <- 

Fm(x)  =  sin  — ^i  8.6 


Hn(y)  =  sin  -2g5  8.7 
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where  a  and  b  are  the  lengths  of  the  edges  along  x  and  y 
directions,  respectively. 

By  substituting  Equations  8.6  and  8.7  into  Equation 
8.5,  one  obtains 


4  4  2  2  4         4  4 

Dll  5LJ_  +  2(d12  +  2D66)  2L^J_  +  D22  £L^  =  xk 

a  a  b  b         8>g 


Letting  b=a  for  a  square  plate  and  considering  only 
the  first  mode  of  vibration  (i.e.,  m=n=l),  one  can 
obtain  X-^   as 


4 

Xl  "  ^  (D11  +  2D12  +  4D66  +  D22>  8'9 

a 


The  natural  circular  frequency  of  the  first  mode,  co^,  is 
defined  as 


ax   =  (-—)*  8.10 


Then  the  first  mode  solution  of  Equation  8.2   is 


r,  .      ttx       .      iry      lcoi  t  0  .  n 

W   =    aii    sin  —   sin  -^-  e     1  8.11 

11  a  a 
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If  the  viscoelastic  behavior  of  the  materials  is 
considered,  the  complex  form  of  Dj^  should  be  included, 
i.e. 

»*  -"2   IS  (DU  +  2DU  +  4D66  +D22)5S  8-12 

a 

After  separating  the  real  and  imaginary  parts  and 
neglecting  the  higher-order  terms  of  the  binomial 
expansion  of  the  quantity  on  the  right-hand  side  in 
Equation  8.12,  the  following  expression  is  obtained. 


w*  =  co^  +   i  w£  8.13 


where 

«i   =  -h       4"    (Dll   +    2D12    +   4D66    +  D22)J5  ■    8'14 


2              D"      +    2D"      +    4D"      +   D" 
„         it                   11              12              66            11  Qic 

oj"    =  — _.  — _ 8.15 

23        P(Dil   +   2Di2   +   4D66    +  D22)iS 


Hence,    for  viscoelastic   material,   Equation   8.11   should  be 
rewritten  as 
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W=    (au   sin  »  sin  SL  eM*)    e-»it  8>lfi 

Equation  8.16    is  similar  to   the  mathematical  model  of 
logarithmic  decrement   [42],    and  the   logarithmic  decrement 
6  can  be  approximated  by 


to. 

2tt 


wT  u" 


6"-      -sj^— *  2tt        1  8.17 


1-     (rrf-) 
Wl 


for  light  damping,  i.e.  when  to?/ti)J  <  <1. 

By  using  Equation  2.22,  the  loss  factor  of  this 
system  is  then  determined  as 


n  =  2 


WI  .  D'il   +   2Dl2  +  4D66  +   D22 

*T    Dil  +  2D12  +  4D66  +  D22  8,:L8 


The  value  of  Amn  depends  on  the  stiffness.   High 
stiffness  reduces  the  deflection  at  resonance;  the 
material  damping  can  further  reduce  the  displacement  at 
resonance.   Gibson  [61]  shows  that  flexural  vibration  for 
a  double  cantilever  beam  under  forced  vibration  at 
resonant  frequency  depends  inversely  on  the  product  of 
material  damping,  area  moment  of  inertia  of  cross  section 
and  Young's  modulus.   It  should  be  noted  that  his  result 
is  based  on  the  equation  of  motion  for  free  vibration  and 
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the  boundary  conditions  of  forced  vibration.   For  the  case 
considered  here,  the  product  of  the  area  moment  of  inertia 
and  Young's  modulus  corresponds  to  the  generalized 
stiffness  D,  which  is  defined  as 


5*  =  D'  +  i  5"  8.19 


D'=  D'  +  2D|2  +  4D£6  +  D"22  8.20 


D-  -  D^  +  2DJ2  +  4D£6  +  D»2  8.21 


Equation  8.18  can  be  rewritten  as 

n  =  DH/D'  8.22 

In  order  to  have  small   resonant  deformation,    a  high 
value  of  the  product  of  stiffness  and  material  damping  is 
required.      But  high  specific  stiffness   is   the  major  pro- 
perty of  fiber  composites   to  be  widely  used  in  space 
vehicle  and  aircraft.      Therefore,    the  mathematical    formu- 
lation of  optimization  on  damping  and  specific  stiffness 
is  to  seek   the  maximum  value  of   the  objective   function 

f-L    (p/d,    s/d)  8.23 
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of  two  variables  p/d  and  s/d,  where  f -^  is  assumed  in  the 

form 


D' 


fx  (p/d,  s/d)  =  T1   B —  +  t2 

5' 

o 


D' 

.    n 

P 

«* 

D' 

o 

•    r'm 

D„      ' 

8.24 


or 


f±    (p/d,  s/d)  =  -El  —1  (T,  +  T2  -4-)  8.25 

p    D '  nm 


The  ranges  for  design  variables  p/d  and  s/d  are 

0  <  p/d  *  0.05  s/d  8.26 

25  i   s/d  £  10000  8.27 

where  T^  and  T2  are  weighting  constants,  nm  is  the  damping 
value  of  epoxy,  p  is  the  density  of  designed  composite,  D' 
and  n  of  design  composite  are  defined  in  equations  8.20 
and  8.22,  respectively,  and  D^  and  po  are  the  correspond- 
ing values  of  D'  and  density  of  continuous  graphite  rein- 
forced epoxy  having  the  same  stacking  sequence  and  fiber 
volume  fraction.   If  high  specific  stiffness  is  important 
in  structures,  T^  could  be  chosen  as  a  high  value. 
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Alternatively,  if  small  resonant  deformation  is  the  major 
consideration,  T2  should  be  a  high  value. 

Case  two  is  to  optimize  the  specific  stiffness  and 
material  damping  of  an  orthotropic  square  plate  clamped  on 
all  four  sides  under  free  vibration. 

The  equation  of  motion  is  just  the  same  as  Equation 
8.2.   The  first  mode  solution  is  assumed,  according 
to  reference  [70] 

W  =  An  F(x)  H(y)  T(t)  8.28 

where 

F(x)  =  8i  cos  X-^x  -  B-l  cosh  X-^x  +  sin  X-^x  -  sinh  X-[_x 

8.29 

H(y)  =  B-l  cos  X-^y  -  B]^  cosh  X-^y   +  sin  X^y  -  sinh  \±y 

8.30 


sin  X,a  -  sinh  X,a 

Bi  =  e — t c-t= —  8-31 

1   -  cos  X.a  +  cosh  X,a 


and  Xi  is  the  constant  to  be  determined. 

The  natural  frequency  of  first  mode  vibration  of  a 
clamped  square  plate  was  given  in  reference  [70]  as 


2  , 

cd-L  =  ^  -    [5.14D1;L  +  1.55(2D22  +  4D66)  +  5.14D22] 

a 

8.32 
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By  a   similar  approach   to   that  prescribed  in  case  one,   one 
will   obtain  the  loss  factor  of  the  system  for  case  two    as 


5.14D"      +    1.55     (2D£2+   4D£g)    +    5.14D^2 
n  =     rmoJJ  +  1.55    (2d{2+  4d^6)    +  5.14d22  8*33 


The  mathematical   formulation  of  optimization  for  this 
case   is   to   seek  the  maximum  value  of  the  objective 
function 

f2     (p/d,    s/d)  8.34 

of  two  variables  p/d  and  s/d,   where  f2  is  defined  in 
Equation   8.37 

The   ranges   of  design  variables   are 

0    <   p/d    <    0.05    s/d  8.35 

25    <    s/d    <    10000  8.36 

where 


D' 


f2     (Pl,    s/d)    =   T-l   — -£—   +   T; 

D' 


D' 

•   n 

P 

D' 

o 

•   rim 

0. 

8.37 


and 
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D'    =   5.14D^   +   1.55{2D'      +   4D£g)    +   5.14D^2 


8.38 


D"    =    5.140^    +    1.55(2DJ2    +    4D^)    +    5.14D^2  8.39 


n   =   D"/D' 


8.40 


CHAPTER  9 
RESULTS  AND  CONCLUSIONS 


9.1   Preliminary  Remarks 

The  damping  analyses  presented  in  Chapters  3  to  6  are 
applicable  to  all  kinds  of  fiber  composites,  including 
continuous  fiber  composites,  discontinuous  fiber  compo- 
sites, symmetrically  or  unsymmetrically  laminated  or  lami- 
nar composites,  randomly  oriented  fiber  composites,  etc. 
The  optimization  analysis  on  damping  and  specific  stiff- 
ness presented  in  Chapter  8  is  based  on  orthotropic 
material.   However,  a  similar  approach  is  applicable  for 
more  general  anisotropic  materials.   Since  most  widely- 
and  practically-used  fiber  composites  are  symmetrically 
laminated  fiber  composites,  the  numerical  analysis  of 
this  study  is  concentrated  on  certain  kinds  of  symmetri- 
cally laminated  fiber  composites,  unidirectional  laminar 
composites,  and  in-plane  randomly  oriented  short-fiber 
composites. 

As  mentioned  before  (see  Chapters  3  and  5),  damping 
nx  is  defined  as  Ex/Ex  (or  Dx/Dx) .   Since  both  Ex  and  E£ 
(or  Dx  and  Dx)  are  functions  of  Ef,  Em,  s/d,  G,  P  (the 
longitudinal  distance  between  fiber  tips),  V^  and  riff  etc., 
the  variations  of  nx  and  Ex  (or  Dx)  may  not  follow  the 
same  pattern.   Consequently,  in  the  numerical  results, 
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three  dependent  variables  (i.e.,  n'  ,  E'  ,  and  E^)  are  dis- 
cussed  in  laminar  composites,  while  more  dependent 
variables  are  involved  in  other  composites.   Among  the 
independent  variables,  most  interest  is  focused  on  Ef,  E  , 
s/d,  and  0.   P  is  treated  as  an  independent  variable  only 
for  optimization,  and  Vf  is  set  as  a  constant  in  most 
cases. 

The  numerical  results  presented  in  this  chapter  are 
based  on  the  following  assumptions: 

1.  The  matrix  epoxy  behaves  as  a  linear,  isotropic, 
viscoelastic  material  [15]. 

2.  All  fibers  are  also  linear  viscoelastic  materials,  and 
damping  of  fibers  is  independent  of  direction. 

3.  The  dependence  of  modulus  on  frequency  is  not 
considered  at  the  present  time. 

9.2   Damping  and  Stiffness  of  Unidirectional 
Fiber  Composites 

The  numerical  results  of  damping  and  stiffness  analy- 
sis presented  in  Chapter  3  are  discussed  in  this  section. 
It  should  be  noted  that  the  modified  Cox's  model 
(i.e.  Equations  3.8  and  3.13)  is  used  to  calculate  the 
longitudinal  modulus  of  composite  materials. 

Figures  9.1,  9.2,  and  9.3  present  the  nondimensional 
plots  of  E^/E^  ,  E^/E"m,  and  nx/nm  of  graphite  epoxy 
composite  as  functions  of  fiber  aspect  ratio  s/d,  using  0 
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(the  angle  between  fiber  direction  and  loading  direction) 
as  parameter.   For  small  angle,  say  0<1O°,  the  storage 
modulus  sharply  increases  as  fiber  aspect  ratio  increases; 
the  loss  modulus  reaches  its  maximum  around  s/d=130  and 
then  reduces  as  s/d  increases,  while  damping  apparently 
decreases  as  s/d  increases.   Comparing  0°  curves  with  the 
30°  curves  in  those  drawings,  one  can  notice  that  when 
0=0°  low  damping  at  large  s/d  is  caused  mainly  by  high 
storage  modulus.   Comparing  the  60°  curves  with  the  90° 
curves  in  those  drawings,  it  is  recognized  that  for  large 
angles,  the  differences  in  loss  moduli  have  more  influence 
on  damping.   Alternatively,  Figures  9.4,  9.5,  and  9.6  show 
the  nondiraensional  plots  of  F^  /E^  ,  E'x  /E'^  ,  and  nx/nm  of 
graphite  epoxy  composites  as  functions  of  0,  using  s/d  as 
parameter.   It  is  clearly  observed  that  the  shorter  the 
fiber,  the  stronger  is  the  dependence  of  Ex/E^  ,  Ex/E^  ' 
and  nx/nm  on  s/d  for  small  angles,  say  0<1O.   When  0>45, 
there  is  not  any  change  of  F^/E^  ,  Ex/E'^  ,  and  nx/nm 
caused  by  changing  s/d.   In  addition,  Figure  9.5  shows 
that  the  maximum  value  of  E^/E^  occurred  between  0=0° 
and  0=15°  for  fiber  aspect  ratio  ranging  from  25  to  10000. 

Similar  properties  of  unidirectional  boron  epoxy 
composite  and  glass  epoxy  composites  are  observed.   But, 
for  0=0,  the  maximum  values  of  Ex/E^  occur  around  s/d=500 
and  s/d=50  for  boron  epoxy  composite  and  glass  epoxy 
composite,  respectively.   Also,  the  damping  of  boron  epoxy 
composite  is,  in  general,  higher  than  that  of  graphite 
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epoxy  composite.    This  is  because  boron  itself  has  higher 
damping  than  graphite  does. 

The  corresponding  drawings  of  Figures  9.1,  9.2,  and  9.3 
for  unidirectional  Kevlar  epoxy  composite  are  given  in 
Figures  9.7,  9.8,  and  9.9,  respectively.   It  is  observed 
that  Kevlar  epoxy  composite  has  the  highest  material  damp- 
ing among  those  four  kinds  of  composite.   This  is  due  to 
the  fact  that  Kevlar  (aramid)  is  a  kind  of  polymer  which 
usually  has  high  damping.   It  is  interesting  to  notice 
that  E"  /E"  of  Kevlar  epoxy  composite  monotonically 

*\  ill 

increases  with  respect  to  s/d  when  0=0°.  All  other  proper- 
ties of  Kevlar  epoxy  composite  are  similar  to  those  of 
graphite  epoxy  composite. 

Figure  9.10  presents  the  plots  of  E^/E^,  E£/E£,  and 
nv/nm  of  Kevlar  epoxy  composite  as  functions  of  0  with 

.X.    Ill 

s/d=100.  This  figure  shows  that  high  damping  of  such 
material  is  induced  by  the  fact  that  the  values  of  E^/E^ 
and  E£  /Ej£  are  very  close.   The  damping  of  Kevlar  epoxy 
composite  is  insensitive  to  angle  change  because  the 
reduction  rates  of  E^/E^  and  E^/E^  are  very  close.   Figure 
9.11  shows  the  plots  of  E^/E^  and  E£/Em  ,  and  nx/nm  of 
Kevlar  epoxy  composite  and  graphite  epoxv  composite  as 
functions  of  0,  while  s/d  is  kept  as  100.   This  suggests 
idea  that  if  a  hybrid  composite  (graphite  epoxy  and  Kevlar 
epoxy  composite)  is  made,  damping  of  this  hybrid  composite 
will  be  higher  than  that  of  graphite  composite,  while  the 
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stiffness  of  the  hybrid  composite  will  be  greater  than  that 
of  Kevlar  epoxy  composite. 

Figure  9.12  presents  E^/E^  of  graphite  epoxy  composite 
as  function  of  0  for  different  fiber  volume  fractions  Vf, 
when  s/d=100.   Differences  in  E^/E^  caused  by  different 
volume  fractions  decreases  when  angle  increases.   It  is 
also  observed  that  an  apparent  reduction  in  E^/E^  occurs 
when  angle  increases  from  0°  to  15°.   Figure  9.13  exhibits 
the  plots  of  nx/nm  for  graphite  epoxy  composite  versus  0 
for  different  fiber  volume  fractions,  when  s/d=100.   The 
peak  value  of  damping  moves  from  about  50°  to  30°  when  the 
volume  fraction  changes  from  0.1  to  0.6. 

Three-dimensional  plots  using  E^  and  nx  as  vertical 
axis,  and  s/d  and  0  as  two  base  axes  of  Kevlar  epoxy 
composites  are  given  in  Figures  9.14  and  9.15, 
respectively.   The  associated  contour  curves  are  presented 
in  Figures  9.16  and  9.17. 

9.3   Damping  and  Stiffness  of  Randomly  Oriented 
Short-Fiber  Composites 

Figure  9.18  gives  the  nondimensional  plots  of  nr/nm/ 

EL/E'  ,  and  E"  /E"mof  glass  epoxy  composite  as  functions  of 

s/d.   It  is  observed  that  damping  reaches  its  maximum  value 

at  small  values  of  s/d  and  decreases  as  s/d  increases.  The 

storage  modulus  E^./E^  monotonically  increases  as  s/d 

increases.   However,  loss  modulus  E£/E^  almost  remains 

unchanged.   In  Figures  9.19  and  9.20,  nondimensional  plots 

of  nr/nm,  E'r/E^  and  E£/E'^  for  graphite  epoxy  composite 
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and  Kevlar  epoxy  composite  are  presented  separately.  The 
variations  of  nr/nm  and  E^/E^  for  Kevlar  epoxy  composite 
or  graphite  epoxy  composite  are  similar  to  those  of  glass 
epoxy  composite.   But,  for  Kevlar  epoxy  composite,  E^'/E^ 
increases  as  s/d  increases,  while  graphite  epoxy 
composite,  E"/E^  has  its  maximum  value  around  s/d=120. 
Figures  9.21,  9.22  and  9.23  show  similar  plots  of  nGr/nGm, 
G^/G^  i    and  G^'/G^  as  functions  of  s/d  for  glass  epoxy 
composite,  graphite  epoxy  composite,  and  Kevlar  epoxy 
composite,  respectively.   Variations  of  riGr/nGm,  G^/G^ 
and  G£/G^  are  almost  the  same  as  those  of  nr/nm,  Er/Em, 
and  E"/E"  of  corresponding  composite,  respectively. 

Figure  9.24  shows  the  nondimensional  plots  of  nr/nm 
as  function  of  fiber  volume  fraction  Vf,  using  s/d  as 
parameter.   It  is  observed  that  nr/nm  decreases  as  Vf 
increases  for  small  s/d  as  well  as  for  large  s/d. 

Three-dimensional  plots,  using  E^.  and  nr  as  vertical 
axis,  and  Ef  and  s/d  as  two  base  axes,  of  four  kinds  of 
fiber  composite  are  given  in  Figures  9.25  in  9.26,  sepa- 
rately.  The  associated  contour  curves  are  presented  in 
Figures  9.27  and  9.28. 

9.4   Damping  and  Stiffness  of  Laminated  Composite — 
Laminated  Plate  Theory  Approach 

The  numerical  results  of  the  damping  and  stiffness 

analysis  described  in  Chapter  5  are  presented  in  this 

section.   Five  types  of  lamination  (quasi-isotropic 

lamination,  angle  ply  lamination,  cross  ply  lamination, 
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unidirectional  lamination,  and  a  special  kind  of  lamination 
[0n/-04rJs)  of  graphite  Epoxy  composites  are  considered  in 
numerical  examples. 

Figure  9.29  exhibits  nondimensional  plots  of  D^  and 
F^ll  (flexural  normal  damping)  for  [02/902/452/-452] s, 
[452/02/-452/902]s,  and  [452/-452/902/02] s  as  functions  of 
s/d.   It  should  be  noted  that  all  those  three  kinds  of 
lamination  have  the  same  values  of  in-plane  longitudinal 
stiffness  A-q  and  in-plane  longitudinal  damping  jTin*   But 
Figure  9.29  shows,  as  expected,  that  for  the  same  s/d, 
[02/902/452/-452]  s  has  the  highest  value  of  D-J^  and  the 
lowest  value  of  Fnn/  [452/-452/902/02] s  has  the  lowest 
value  of  Diq  and  the  highest  value  of  pin/  while  [452/02/ 
-452/902]   has  almost  the  average  values  of  D^  and  of 
Fnll  °^  ^e   otner  two  kinds  of  lamination.   Alternatively, 
Figure  9.30  shows  the  nondimensional  plots  of  Dgg  and  pilgg 
(twisting  damping)  for  those  three  kinds  of  lamination  as 
functions  of  s/d.   In  this  drawing  for  same  s/d,  [02/902/ 
452/-452]s  has  the  lowest  value  of  Dgg  and  the  highest 
value  of  pHgg/  [452/-452/02/902] s  has  the  highest  value  of 
Dgg  and  the  lowest  value  of  Fngg/  and  [452/02/-452/902] s 
still  remains  around  the  average  values  of  those  other  two 
kinds  of  lamination.   Figures  9.29  and  9.30  actually  show 
us  how  to  design  the  quasi-isotropic  laminated  plate  for 
either  maximum  damping  or  maximum  stiffness  purpose. 

Figures  9.31  and  9.32  show  the  nondimensional  plots  of 
Ej^  and  jnn  (in-plane  longitudinal  damping)  of  angle  ply 
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laminated  graphite  epoxy  composite  as  function  of  G,  using 
s/d  as  a  parameter.   E^  increases  with  s/d  for  small 
angles,  say  0<2O°.   The  variation  of  Ej^  with  s/d  reduces 
as  G  increases,  and  very  small  differences  exist  for  0>6O°. 
The  jln  distinctively  decreases  as  s/d  increases  for 
small  angles  (say,  G<30°),  while  the  variation  of  j^n 
with  s/d  reduces  as  G  increases.   Figures  9.31  and  9.32 
are  comparable  to  Figures  9.4  and  9.6  correspondingly.   The 
reduction  in  jn^  (or  the  increment  in  Eq)  of  angle  ply 
laminated  composite  is  not  as  sharp  as  that  of  unidirec- 
tional composite,  when  s/d  is  increased  and  the  angle  is 
small . 

In-plane  shear  stiffness  Egg  and  in-plane  shear 
damping  jHgg  of  angle  ply  graphite  epoxy  composite  as 
function  of  G,  using  s/d  as  parameter  are  shown  in  Figures 
9.33  and  9.34,  separately.   Maximum  Egg  and  minimum  rTigg 
occur  at  G=45°.   Larger  dependences  of  Egg  and  jHgg  on  G 
are  observed  for  larger  s/d.   The  reason  why  D'i  i,  /D'  , 

D66/DGm'  Fnll/rW  and  Fn66/nGm  are  not  Presented  in 
drawings  is  that  those  values  are  very  close  to  the  values 

of  Eil/Em  '  E66/Gm  '  I^ll^m'  and  i^e^Cm  separately. 

Figures  9.35-9.38  present  the  nondimensional  plots  of 

Dll'  F^ll'  D66'  and  Fn66  °f  graphite  epoxy  composite 
laminated  with  four  kinds  of  lamination  ([45/-45]4  , 
[452/02/-452/902]s,  [0/90]4gf  and  [0g]s)  as  functions  of 
s/d.   It  is  observed  that  for  the  same  s/d,  [45/-45]4s  has 
the  highest  values  of  Dgg  and  Fnii  and  the  lowest  values 
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of  D-^  and  Fn66.   The  values  of  D£g  and  Fn66  for  [0/90]4gf 
and  for  [0g]s  are  same,  and  those  values  are  almost 
independent  of  s/d.   This  is  because  ET  and  GLT  of 
fiber  composites  are  assumed  to  be  indendpent  of  the  fiber 
length.   D]^  of  [0g]g  is,  as  expected,  higher  than  D-^  of 
[0/90]4g  at  the  same  s/d,  but  Fnn  of  those  two  kinds  of 
lamination  are  surprisingly  close.   The  D^  ,  pniii  D^6 

and  prigs  of  ^ 452/°2/~452//902^  4s  are  Just  somewhere  near 
the  middle  range  of  those  of  [45/-45]4s  and  [0/90] 4  . 

The  characteristic  property  of  [©n/-04n]s  is  that 
its  D[6  and  D^6  are  very  close  to  zero.   For  example, 
D16//Dil  or  D26  /D'n  of  this  lamination  for  graphite  epoxy 
composite  is  less  than  0.02  for  all  angles,  while  those  of 
angle  ply  lamination  are  less  than  about  0.11  for  all 
angles.   Numerical  results  show  that  Dj^,  f^II'  d66'  an(* 
Fn66  of  f02//-G8-'s  are  verY  close  (difference  within  2%)  to 
those  of  [0/-G]5s.   But  Ej^  and  E£6  (or  1r\11   and  xn66) 
are  in  general  less  (or  greater)  than  those  of  [G/-G]r 
(maximum  differences  ranging  from  13%  to  30%) . 

9.5   Damping  and  Stiffness  of  Laminated  Fiber  Composite — 

Energy  Approach 

There  are  twenty-one  different  kinds  of  material 

damping  induced  from  Cft  /C{j  (i,  j  =  1,  2  .  .  .  6)  for 

anisotropic  materials.   In  the  energy  approach,  material 

damping  of  structure  is  just  the  ratio  of  the  sum  of 

h   efj  CVj  to  the  sum  of  \    efj  Cj_j.   Therefore,  different 

values  of  material  damping  could  be  observed  for  one 
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material  under  different  kinds  of  loading  and  boundary 
conditions.   In  the  author's  opinion,  damping  analysis 
through  the  energy  approach  is  more  acceptable  than  lami- 
nated plate  approach,  because  the  energy  approach  consi- 
ders the  influence  of  all  stresses  on  the  dissipated 
energy,  while  laminated  plate  approach  considers  one 
stress  at  a  time.   The  results  of  damping  analysis  through 
the  energy  approach  presented  in  Chapter  6  are  discussed 
in  this  section. 

Figure  9.39  shows  the  difference  between  material 
damping  through  the  energy  approach  and  through  the  lami- 
nated plate  approach  for  a  60x60x10  cross  plied 
(O5/9O5/O5)  graphite  epoxy  composite  under  uniform  in- 
plane  tension  load  along  x  direction.   The  difference  is 
mainly  caused  by  the  fact  that  there  exists  some  stress 
other  than  ax;  for  example,  txz  is  not  negligible  in 
certain  regions.   The  difference  increases  when  the  value 
of  s/p  is  increased.   In  other  words,  continuous  fiber 
composites  have  more  apparent  difference  of  damping  values 
through  those  two  approaches  than  discontinuous  (short- 
fiber)  composites  do. 

Figure  9.40  presents  the  influence  of  plate  thick- 
ness and  of  the  lamination  on  the  damping  value  by  the 
energy  approach  under  in-plane  uniform  tension  load.   It  is 
noticed  that  either  increasing  the  number  of  interfaces  or 
reducing  the  total  plate  thickness  will,  in  general,  reduce 
the  difference  between  damping  values  of  those  two 
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approaches.   It  should  also  be  noted  that,  for  an  isotro- 
pic material  under  tensional  load  along  x  direction,  the 
damping  value  of  the  energy  approach  is  same  as  the  value 
of  tH^i  through  the  laminated  plate  theory  approach. 

When  a  flexural  bending  load  is  considered,  a  pair  of 
couple  forces  are  applied  around  one  edge;  while  a  fixed 
boundary  condition  is  presumed  at  the  opposite  edge. 
Figure  9.41  shows  that  reducing  the  total  thickness  makes 
the  damping  values  of  the  energy  approach  closer  to  those 
of  the  laminated  plate  approach.   It  is  also  indicated  by 
Figure  9.41  that  for  the  flexure  bending  load  case, 
increasing  the  total  number  of  interfaces  would  somehow 
help  those  two  approaches  to  obtain  more  consistent 
damping  results. 

9.6   Experimental  Results  of  Damping 
and  Stiffness 

Figure  9.42  presents  some  experimental  results  by  the 
impulse  hammer  technique  [35,  36]  about  material  damping 
and  stiffness  of  graphite  epoxy  composite  beams.   Experi- 
mental results  of  stiffness  are  very  close  to  those  of  the 
analytical  results,  while  experimental  results  of  damping 
are  not  so  close  to  the  analytical  results.   A  similar 
trend  is  also  presented  in  Figure  9.42  where  experimental 
results  [35,  36]  about  damping  and  stiffness  of  off-axis 
unidirectional  continuous  graphite  epoxy  composite  beams 
are  compared  with  analytical  results.   Table  9.1  shows  the 
experimental  results  of  damping  and  of  the  first  mode 
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nature  frequency  (this  quantity  is  proportional  to  the 
stiffness  of  the  structure)  of  five  cross  ply  short-fiber 
glass  epoxy  composite  plates.   It  is  observed  that  high 
value  of  s/p  (either  increase  fiber  length,  s,  or  reduce 
the  gap  between  fiber  tips,  p)  will  increase  the  nature 
frequency  and  decrease  the  damping. 

9.7   Optimization  of  Damping  and  Specific  Stiffness 

of  Fiber  Composites 

According  to  the  prescription  about  optimizing  damp- 
ing in  Chapter  8,  some  optimum  designs  are  discussed  in 
this  section  for  graphite  epoxy  composites  and  for  hybrid 
composites  keeping  Vf  equal  to  0.65.   The  design  vari- 
ables include  the  fiber  aspect  ratio,  s/d,  and  the  ratio 
of  the  longitudinal  distance  between  fiber  tips,  p,  to  the 
fiber  diameter,  d. 

Table  9.2  shows  the  influence  of  weighting  constants, 
T-j_  and  T2,  on  the  optimum  design  for  0/90/0  graphite  epoxy 
composites  under  situations  of  case  one.   By  increasing 
T2,  s/p  is  reduced;  while  by  increasing  T-,,  s/p  is 
increased.  These  imply  that  to  have  smaller  resonant 
deformation,  larger  distance  between  fiber  tips  is 
preferred;  while  to  have  higher  specific  stiffness, 
smaller  distance  between  fiber  tips  is  required. 

Table  9.3  shows  the  optimum  design  for  cross  ply 
graphite  epoxy  composites  and  hybrid  composites  under  case 
one  situation  and  having  the  same  total  thickness,  with 
Tj^l  and  T2=5.   There  is  no  difference  in  the  optimum 
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values  of  design  variables  s/d  and  p/d  for  0/90/0  and 
0/90/0/90/0  graphite  epoxy  composites.   For  the  same  total 
thickness,  with  1^=1  and  T2=5,  discontinuous  fibers  give  a 
high  value  of  objective  function  in  the  stacking  sequence 
0G/90G/0K/90K/s;  while  continuous  fibers  give  high  value 
of  the  objective  function  in  the  stacking  sequence 
0K/90K/0G/90G/90G/s  where  the  superscripts  G  and  K  repre- 
sent graphite  and  Kevlar,  respectively,  and  the  fiber 
composites  having  s/p=2593.  can  be  treated  as  continuous 
fiber  reinforced  composites.   When  comparing  the  value  of 
objective  function  of  each  optimum  design  for  each  case, 
one  can  discover  that  0K/90K/0G/90G/s  has  the  highest 
value  of  objective  function.   This  indicates  that  the  con- 
tinuous hybrid  fiber  composite  plate  with  the  stacking 
sequence  0K/90K/0G/90G/s  is  the  optimum  design  for  case 
one. 

Figure  9.44  shows  the  contour  curves  of  objective 
function  of  case  one  vs  s/d  and  s/p  (keeping  Ti=l,  T2=5). 
It  is  observed  that  there  probably  exists  an  acceptable 
region  as  optimum  design.   In  fact,  the  term  (s/p)ave  in 
Tables  9.2,  9.3  and  9.4  is  the  average  value  of  ratio  s  to 
p  at  all  points  whose  objective  value  is  within  99.95%  and 
100.05%  times  of  the  objective  value  of  the  optimum  design 
shown  in  s/d  and  p/d  items.   Numerical  results  show  this 
region  ranging  from  s/d=2000  to  s/d  over  9000  for  case 
one.   That  means,  as  long  as  s/p  is  kept  equal  to  53.1  and 
s/d  over  2000,  one  will  get  an  acceptable  optimum  design 
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(when  Tx=l  and  T2=5)  for  0/90/0  or  0/90/0/90/0  cross  ply 
graphite  epoxy  composite  plates  under  case  one  conditions. 

The  optimum  designs  of  case  two  (see  Table  9.3)  is 
very  close  to  those  of  case  one,  although  D  for  those  two 
cases  are  different.   It  is  because  D-^  and  Dgg  are  too 
small  relative  to  D1;L  and  D22  (see  Equation  8.20  and 
Equation  8.38) . 

9.8   Conclusions 
Since  the  material  damping  and  stiffness  of  continuous 
or  discontinuous  fiber  composites  depends  upon  many  quanti- 
ties such  as  stiffness  ratio  Ef/Em,  fiber  aspect  ratio 
s/d,  longitudinal  distance  beetween  fiber  tips  p,  loading 
(or  fiber)  angle  0,  fiber  volume  fraction  Vf,  fiber  and 
matrix  damping  values,  and  stacking  sequence,  it  is  not 
feasible  to  discuss  the  material  damping  and  stiffness  by 
considering  all  the  possible  variations  of  all  the 
independent  quantities  simultaneously.   However,  based 
upon  the  numerical  results  presented  in  this  study,  we  can 
make  the  following  conclusions: 
A.   About  Unidirectional  Fiber  Composites 

1.  High  damping  could  be  achieved  by  small  0,  small 
s/d,  and  high  p;  while  high  stiffness  could  be 
obtained  through  high  Vf,  small  0,  large 

s/d,  and  low  p. 

2.  Damping  of  discontinuous  fiber  composites  is  no 

less  than  that  of  continuous  fiber  composites. 
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3.   For  small  angle,  high  s/d  gives  high  stiffness; 
while  low  s/d  brings  high  damping.   For  large 
angles,  damping  and  stiffness  are  independent  of 
s/d. 

B.  About  In-Plane  Randomly  Oriented  Fiber  Composites 

1.  Relative  to  unidirectional  fiber  composites,  in- 
plane  randomly  oriented  fiber  composites  have 
lower  stiffness  at  the  price  of  keeping  high 
damping  in  all  directions  parallel  to  the  plane. 

2.  Variations  of  shear  damping  and  shear  modulus  vs 

s/d  are  very  close  to  those  of  normal  damping 
and  normal  stiffness,  respectively. 

3.  Influence  of  Vf  on  normal  damping  becomes  less 
when  s/d  is  smaller. 

C.  About  Laminated  Fiber  Composites 

1.  For  different  kinds  of  lamination,  high  stiffness 
lamination  has  low  damping  and  vice  versa. 

2.  The  properties  of  damping  and  stiffness  of  angle 
ply  lamination  are  very  similar  to  those  of 
unidirectional  fiber  composites,  when  0  is  the 
same. 

3.  Damping  of  cross  ply  lamination  is  close  to  that 
of  aligned  fiber  composite  (0=0°).   But  the 
cross  ply  keeps  equal  stiffness  in  both  x  and  y 
directions  rather  than  high  stiffness  in  x 
direction  and  low  stiffness  in  y  direction  of  the 
aligned  fiber  composite. 
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D.    About  the  Energy  Approach  and  Laminated  Plate  Theory 
Approach 

1.  Damping  analysis  through  the  energy  approach  is 
loading  condition,  boundary  condition  and 
geometry  (of  the  body)  dependent.   It  gives  more 
reliable  results  on  damping,  while  the  laminated 
plate  theory  approach  can  give  quick  estimation 
about  damping. 

2.  Fortunately,  the  variations  of  the  material 
damping  vs  those  variables  in  these  two 
approaches,  in  general,  follow  a  similar  pattern, 
if  the  dimension  in  the  thickness  direction  is 
one  order  less  than  the  dimensions  in  the  other 
two  directions. 

E.    About  Optimization 

Optimumization  of  damping  and  specific  stiffness 
depends  on  the  preference  of  high  specific  stiffness 
or  high  damping  (T^  or  T2),  material  property  of 
fiber,  the  stacking  sequence  for  hybrid  composite, 
and  boundary  conditions.   A  continuous  fiber  rein- 
forced hybrid  orthotropic  composite  plate  with  stack- 
ing sequence  0K/90K/0G/90G  is  the  optimum  design  for 
clamped  (or  simply  supported)  square  plate  under  the 
free  vibration. 
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Figure  9.14:   Three-Dinensional  riots  of  E'/E'  vs  9  and  s/d 
for  Unidirectional  Kevlar  Epoxy  Composites 
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Figure  9.16:   Contour  Curves  of  E'/E'  vs  0  and  s/d  for 
Unidirectional  Kevlar  Epoxy  Composites 
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Figure  9.25:   Three- Dimensional  Plots  of  E'/E'  vs  s/d  and 
El   for  Randomly  Oriented  Fiber  Composites 
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Figure  9.26:   Three-Dimensional  Plots  of  n  /n  vs  s/d  and  EL 
for  Randomly  Oriented  Fiber  Composites 
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Contour  Curves  of  E^/E^  vs  s/d  and  E£  for 
Randomly  Oriented  Fiber  Composites 
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Figure  9.44: 


Contour  Curves  of  Objective  Function  for  Case 
One  vs  s/d  and  s/p  for  Graphite  Epoxy 
Composite  Plate 
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Table  9.1:   Experimental  Results  of  05/905/0c  Glass 
Epoxy  Composite  Plates 


s  (mm) 

p  (mm) 

d  (mm) 

Vf 

fr(Hz) 

n 

19. 

1.59 

.009 

.558 

92.8 

.0042 

25.4 

1.59 

.009 

.565 

95.3 

.0041 

38.1 

1.59 

.009 

.575 

97.1 

.0040 

25.4 

0.05 

.009 

.579 

109.8 

.0036 

Table  9.2:   Influence  of  Weighting  Constants  on  Optimum 
Design  for  Case  One 


Lamination   Tl     T2    s/d     p/d    (s/p)ave    g,     j] 

pDl 

0/90/0 1     3    9022.    121.    74.4 1.27    .47 


0/90/ 

1 

5 

9022. 

170. 

53.1 

1.19 

.52 

0/90/0 

1 

10 

9008. 

214. 

42.1 

1.12 

.56 

0/90/0 

0 

5 

9020. 

267. 

33.8 

1.06 

.60 
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Table  9.3: 


Optimum  Designs  of  Cross  Ply  Composite 
Plates  for  Case  One 


~s7H    p7B    Ci7pT 


Lamination 


ave 


D'p* 
pDi 


H 

n 


m 


(0/90/0) 


(0/90/0/90/0) 


5   9022 


5   8994 


170 


169 


53.1 


53.1 


1.19 


1.19 


0G/90G/0K/90K/s   1    5 9007.    152.    59.4      1.29 

0K/90K/0G/90G/s   1    5   9993.    3.2    2593.1     1.11 


52 

_52 

53 

69 


Table  9.4:   Optimum  Designs  of  Cross  Ply  Composite 
Plates  for  Case  Two 


Lamination 


Tl   T2 


T/d        pTH  TsTpT 


ave   D'p. 


pDS 


0K/90K/0G/90G/s    1 


9993. 


3.2 


2747.8   4.80 


m 


(0/90/0)G 

1 

5 

9010. 

170. 

53.0 

5.10 

.46 

(0/90/0/90/0)° 

1 

5 

9020. 

170. 

53.0 

5.10 

.46 

0G/90G/0K/90K/s 

1 

5 

9010. 

152. 

59.3 

5.53 

.48 

.65 


APPENDIX  A 
HYSTERESIS  LOOP  OF  VISCOELASTIC  MATERIALS 


This  appendix  shows  that  hysteresis  loop  could  also  be 
observed  in  a  viscoelastic  material. 

For  a  viscoelastic  body  under  a  cyclic  load  a  at 
circular  frequency  a), 


a   =  oQ   sin  wt  A.l 


strain  e  of  this  body  lags  a  phase  angle  $  behind  tress, 
i.e. , 


e  =  eQ  sin(tot  -  $)  A. 2 


where  cQ   and  aQ   are  the  magnitudes  of  strain  and  stress 
respectively,  and  e0,  in  general,  is  proportional  to  aQ. 


o0  =  E  e0  A. 3 


where  E  is  the  elastic  Young's  modulus.   Equation  A.l  can 
be  rewritten  as 


a  =   E  eQ  sin(wt  -  $  +  $)  A. 4 
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that   is 


a   =   E    e      cos$    sin(wt   -    <J>)    +   E    gq   sin<i>    cos(ut   -    $) 

A. 5 


By  the  definitions  of  storage  and  loss  moduli,  Equation  A.5 
induces  to 


a  =  E'  eQ  sin(wt  -  *)  +  E"  eQ  cos (wt  -  $) 


From  Equations  A. 6  and  A. 2 


A. 6 


sin(tot   -    $)    =  — - —  A. 7 

eo 


COS(wt    -    <J>)     =    ^snrz A-8 

E   eo 


Using  Equations  A. 7  and  A. 8,  one  will  have 


2       2 

(  gFnrE'£  )   +  (— )   =1  A. 9 

E  eo         eo 


Fiaure  A.l  is  the  sketch  of  Equation  A. 9 
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Figure  A.l:    Hysteresis  Loop  of  a  Viscoelastic  Material 

Equation  A. 9  shows  that  hysteresis  loop  is  also  presented 
in  viscoelastic  material. 

The  dissipated  energy  UD  during  a  cycle  is  the  area 
within  the  ellipse. 


(Un)     =  n  £n  E"  e, 
O   eye      o     < 


A. 10 


The  maximum  energy  stored  in  the  body  Ug  is  the  area  of 
the  triangle  OAB 


Therefore, 
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Us   '   5    eo   E'    eo 


A. 11 


E"  (Up)cyc 

TT  "      2ttU„ 


A. 12 


APPENDIX  B 
EXPRESSIONS  OF  6's 


This  appendix  presents  the  expressions  of  8/  Bs/2, 
and  8*s/2.   The  8  used  in  Section  3.1  is  originally 
defined  in  [20]  as 

6  =  [-S5-   --i )h  B.l 

Ef       TJ  Ln(R/Yo) 

where  Gm  is  the  shear  modulus  of  the  matrix,   Ef  is  the 
Young's    modulus    of    the   fiber,    y0   is    tne    radius   of   fiber, 
and  R  is   the  outside   radius   of   the  matrix,    which   surrounds 
the   fiber   in  a   cylindrical    shape,    as   in  Figure  3.2   when 
0=0°.      The   ratio  R/y0  may  be  expressed  in  terms  of  the 
fiber  volume    fraction  Vf    [31].      For  example,    for  a 
hexagonal   array 


R      -,  ! B.2 


y 


o        2      3V 


f 


and  for  a  square  array, 
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Y 


B.3 


o    4  V 


f 


Consequently,  for  a  square  array, 


£--5   -^—  [ 1 lk  B.4 

f     Ln(WJ} 


where  s  and  d  are  the  length  and  diameter  of  the  fiber, 
correspondingly. 

When  a  viscoelastic  material  is  considered, 


13  *  s  _  s    r4     ,*5      rm         Gm   ,  ^ 

1 5     Ln(  !  )  Ef  (1  +  i  nf} 

V  B.5 


where  nGm  =  G£  /G^  B.6 

Gm  '  G'man<^  If  are  defined  in  Equations  3.26  and  3.29. 
Since  the  loss  factors  are  very  small,  high  order  terms 
such  as  nf  and  n^  nGm  could  be  neglected.   Using  a 
binomial  expansion  reduces  Equation  B.5  to 


^£  =  6s  u  +  .  l  (ncm  _  nf)]  B>7 
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Similarly,  by  Taylor's  series  expansion  and  neglect  of 


higher  order  terms,  one  gets 


^„u  6*s      .  8s  ,  .  6s  nGm  "  T]f  1  _  B 

tanh  -j-  -  tanh  j-  +   i  2 2 _-        B.8 

coshw— 


APPENDIX  C 
AN  ALTERNATIVE  WAY  OF  DETERMINING  Er  AND  Gr 


Another  derivation  of  Er  and  Gr  for  in-plane  randomly 
oriented  short-fiber  composite  is  presented  in  this 
appendix. 

Considering  an  in-plane  randomly  oriented  short-fiber 
composite  plate.  Assuming  that  ax  is  the  only  nonzero 
applied  stress,  one  has 


Q,.e   +  Q-,  ~  e   +  Qi  <r  Y    =   o 

11  x  v12   y    *16  'xy     x 

Q,  0  e   +  Q_„  e   +  Q_,  y    =0 

12  x  *22   y    *26  'xy 


C.l 


■16     x 


26   y    66  'xy 


where  Qjh  (1,  k=l,  2,  6)  are  defined  in  Equations  4.1-4.5. 
By  Cramer's  rule 


ex  ■ 


0 
X 

Q12 

°16 

0 

Q22 

52  6 

0 

S26 

he 

Q22   Q26 
°26   566 


C.2 


where 
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A  = 


;11 


;12 


'16 


Q12  Q16 
Q22  C26 
S26   566 


C.3 


According  to  the  definintion  of  Young's  modulus, 


ER  =  T, 


C.4 


Therefore, 


ER  = 


Q22  Q26 
Q26  §66 


C.5 


After  simplifying, 


~2    ~2 
Qll  "  Q12 


JR 


Q 


11 


C.6 


By  Equation  4.6,  Equation  C.6  can  be  rewritten  as 


66 


ER  -  4  Q66  (1  -  —Z ) 


Q 


11 


C.7 


which  is  same  as  Equation  4.11 

A  similar  procedure  by  letting  xxy  be  the  only  non- 
zero applied  stress  can  lead  to  the  same  expression  of 
Equation  4.8  for  Gr. 
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It  should  be  noted  that  the  expression  of  the  denomi- 
nator in  Equation  6.5  is  a  way  to  determine  the  S±±.      Where 
S-,1  is  an  element  of  [S]  ,  and  [S]  is  the  inverse  matrix  of 
[Ql. 


APPENDIX  D 
INVERSION  OF  A  SYMMETRIC  COMPLEX  MATRIX 


This  appendix  presents  a  way  to  get  the  inverse  matrix 
of  a  symmetric  complex  matrix. 

Considering  a  3x3  complex  matrix, 


Al 

+ 

i 

A" 

A  ' 

A3 

+  i  A" 

A  * 
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A" 
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A" 

A6 
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A" 

A5 
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+ 
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where 


*   *   *    *    *    * 
tAl  A2  A3   A4   A5   A6] 


*    *    *    *    *    * 
[A11  A22  A12  A66  A16  A26J 


D.2 


Let  [a  ]  be  the  inverse  matrix  of  [A  ] . 


[a*]  = 


'aj  +  i  a'l 
a3  +  i  a3 


a'  +  i  a 


a'   +  i  a3" 


a2   +  i  a' 


+  i  a 


■   +  i  «5 

;   +  i  a" 


+  i  a 


4 
D.3 


where 
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[a*  a*  a*     a*   a*   a*]  = 


Therefore, 
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[all  a22  a12  a66  a16  a26] 


D.4 


[A*]  [a*]  =  [I] 


D.5 


where  [I]  is  a  unit  matrix: 


Zij  " 


i  =  D 


i  ¥   J 


D.6 


According  to  matrix  operations, 


Iii  =  1  +  i  o  =  (A,'  +  i  A")  (ai  +  i  a")  + 


■11 


1'  v  1 
(Aj  +  AJ) (a'  +  i  aj)  +  (A£  +  i  A£) (a£  +  aj) 


Separating  the  real  and  the  imaginary  parts  reduces 
Equation  D.7  to  the  following  two  equations: 


D.7 


1  =  AjaJ  -  A»a»  +  Ajaj  -  A»a»  +  A- a'-  A»a« 


D.8 


0  =  AjaJ  ♦  A«aJ  +  A'a"  +  A''a'  +  A'a"  +  A»a' 


D.9 
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By  applying  similar  procedures  to  ±12/  1i^'    *22'  J23'  and 
133,  a  system  of  twelve  equations  containing  twelve 
unknowns  (a|  ,  a£  ...  a^  ,  a£ )  is  obtained.   The  solution 
to  this  system  of  equations  determines  [a  ] . 


APPENDIX  E 
ENERGY  EXPRESSION  OF  DAMPING  FOR  LAMINATED  COMPOSITES 


Equation  2.22  indicates  that  the  material  damping,  r\t 
equals  to  the  ratio  of  the  energy  dissipated  in  a  cycle, 
(UD)C   ,  the  product  of  2f,  and  the  maximum  strain  energy 
stored  in  the  body,  Ug,  under  a  periodic  vibration.   Since 
the  maximum  strain  energy  of  a  beam  under  static  load  and 
under  first  mode  vibration  are  equal,  if  the  maximum 
deflection  are  the  same.   If  the  material  damping  of  this 
beam  is  known,  the  energy  dissipated  in  the  beam  within  a 
cycle  of  the  first  mode  vibration  can  be  determined  by 


<Vcyc  =  2u  n  Us  E.l 


where  Ug  can  be  determined  by  calculating  the  strain 
energy  of  the  statiscally-defcrmed  beam  having  the  same 
amplitude  of  the  deflection  as  that  of  the  first  mode 
vibration.    As  it  is  well  known  that  the  dominant  stress 
in  a  beam  under  flexure  vibration  is  normal  stress,  say 
ax.  Therefore,  Us  and  (UD)  cyc  can  be  obtained 
approximately  as 


Us  =  \\   °x  ex  dv  E-2 
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and 


(Vcyc  =  2°  HX   i  ax  gx  dV  E.3 


Where  ev  is  the  normal  strain,  and  nv  is  determined  by 
E"/E'  .   To  be  more  accurate,  all  stresses  and  strains  in 
the  beam  should  be  included  when  calculating  the  maximum 
strain  energy  and  the  dissipated  energy. 

Ungar  and  Kerwin  [56]  derived  the  expressions  for  the 
loss  factor  of  any  series-parallel  array  of  m  viscoelastic 
springs 


m 

.1      (uDj>cyc 
2zl E.4 


n  = 

in 

y      o  tt  r, 


m 
j=l 


1         2 
U  •  =  -  k-   x-  E.5 

SD    2   J    3 


where  Ue^  denotes  the  maximum  strain  energy  stored  in  a 
s  J 

deflection  x^  in  the  jth  lossless  spring  of  stiffness  kj, 
and  (UDj)cyc  denotes  the  energy  dissipated  in  the  jth 
viscoelastic  spring  during  a  cycle.   Since  there  are  six 
different  components  in  the  stiffness  matrix  for  an 
orthotropic  material,  one  can  treat  an  orthotropic 
material  as  a  structure  composed  of  six  springs. 
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Consequently,  the  material  damping  of  a  laminated  n  plies 
composite,  n,  in  the  first  mode  vibration  can  be 
determined  by 


n 

Z  (qUD^  eye 
=   q=1 E.6 


n  =    n 


q=l  4 


where 


?s       I  «j   ej  dv=       j  £k  %   £j  dv 


E.7 


C^j  +  i  Ckj   Hkj  =  C^j  +  i  C£j   =  C*j  E.9 


In  which  V   is  the  volume  of  the  qth  layer  in  the 

Si 

laminated  composite,  C^j  ,  C£j   and  cjj  are  the  storage 
moduli,  the  loss  moduli,  and  the  complex  moduli  for  a 
viscoelastic  material,  respectively. 


APPENDIX  F 
FORMULATION  OF  FINITE  ELEMENT  METHOD 


This  appendix  briefly  prescribes  the  formulation  of 
the  finite  element  method  utilized  in  energy  approach  to 
determine  the  displacement  and  the  strain  fields  for  an 
elastic  body.   The  material  constants  are  also  presented 
here. 

F.l   Stiffness  Matrix 
The  nonzero  elements  of  stiffness  matrix  [C]  for  an 
orthotropic  material  when  its  coordinates  coincide  with 
principal  coordinates  are  given  in  Reference  [57]  as: 

C1;L  =  (1  -  v23  v32)E1  A  F.l 

C12  =  (v12  +  v32  v13)E2  A                        F.2 

C13  =  (v13  +  v12  v23)E3  A                        F.3 

C22  =  (1  -  v13  v31)E2  A  F.4 

C23  =  (v23  +  v21  v13)E3  A                         F.5 

C33  =  d  "  V12  V21*E3  A  F*6 
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C44    -    G23 


F.7 


C55    ~   G31 


F.8 


C66    _   G12 


F.9 


where 


A   = 


1   "    v12v21 


_    v23v32    "    v31v13    "    2VonV,oV1 


'21v32v13 


F.10 


Since  the  material  properties  of  a  laminar  fiber  composite 
with  the  coordinates  parallel  to  and  transverse  to  the 
fiber,  as  Figure  E.l,  can  be  treated  as  transversely 
(in  2,    3  directions)  isotropic, 


Figure  P.l:   Principal  coordinates  of  fiber  composite 
material 


one  has 


El   =  EL 


E2   ~~  ET 


E3   —  ET 


G12  "  GLT 


G13  "  GLT 


■  ET 

'23  "  2  (1  +  vTT.) 


and 


v12  ■  VLT 
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F.ll 


v21 

= 

VLT  ET 

EL 

v23 

= 

Vmm  1 

v32 

= 

Vmm  1 

v13 

- 

VLT 

v31 

= 

VTL 

F.12 
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where  EL,  ET,  GLT,  vLT  are  defined  by  Equations  3.13, 
3.15-3.17,  and  vTT.  is  assumed  to  be  the  Poisson's  ratio 
of  the  matrix.   By  substituting  Equations  F.ll  and  F.12 
into  Equations  F.1-F.9,  the  [C]  matrix  of  transversely 
isotropic  mate-  rial  is  then  obtained.   Equations  F.ll  and 
F.12  show  that  there  are  only  five  independent  material 
constants,  namely,  EL,  ET,  GLT,  vLT,  and  vTTi. 

F.2  Formulations  of  Finite  Element  Method 
A  three-dimensional  elastic  body  bounded  by  surface  S 
is  loaded  by  concentrated  forces  p  on  nodes  and  surface 
traction  t  on  surface  SI.   Assume  that  the  body  force  is 
negligible  and  that  the  whole  system  is  adiabatic  and 
conservative.   According  to  the  principle  of  stationary 
potential  energy  [58],  the  equilibrium  is  ensured  if  the 
total  potential  energy  is  stationary  for  variation  of 
admissible  displacement,  i.e. 

6  (U  -  We)  =  0  F.13 

where  U  is  the  strain  energy,  and  we  is  the  work  done  by 
external  loads.  In  fact,  in  the  elastic  case,  potential 
energy  is  not  only  stationary  but  is  a  minimum  [59]. 

If  there  is  not  any  initial  stress  or  strain,  the 
strain  energy  U  and  the  work  done  by  the  external  load  are 


D  =  \        |   {e}    {o}    dV  F.14 
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and 


m      T 

W~  =  I         {p}   {u}   + 


i=j 


l     l 


ft}   {u}  dA 


F.15 


where 


[e\  [ex  ey  ez  YyZ  YXz  ^xy-' 


F.16 


{a}   =  [ax  a  oz  x   txz  t   ] 


F.17 


T 

|pi  =  rp   p  .  p  •  ] 

\f;      lrxl  cyi    rziJ 

i 


F.18 


{Ui}   =  [uxi  u  i  uzi] 


F.19 


{u}   =  [ux  u  uz] 


F.20 


ft}    =  [tx  ty  tz: 


F.21 


in  which  uv ;   is  the  nodal  displacement  in  x  direction  at 
the  ith  node,  and  PV1-  is  the  concentrated  force  at  the  ith 
node  in  x  direction,  m   is  the  total  number  of  nodes  where 
concentrated  load  is  applied. 
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By  introducing  stiffness  matrix  [C]  and  the  shape 
function  matrix  [N] ,  i.e. 

{o}    =  [C]  {e}  F.22 

{u}  =  [N]  {q}  F.23 

strain  {e}    can  be  expressed  as  function  of  nodal 
displacement  {q}  by  using  strain-displacement  relation 

{e}  =  [L]  {q}  F.24 

Therefore,  Equations  F.14  and  F.15  can  be  rewritten  as 


n  T     T 

U  =  i      {q}   [L]   tC]  [L]  {q}  dV  .  F.25 


m     T  T 

Wp  =   I   {P}   [N]  {q}   +     {t}   [N]   {q}  dA 


e 

i=3     i        1 

sl 


F.26 


By  using  Equations  F.25  and  F.26,  a  system  of  3xn  equa- 
tions is  obtained  from  Equation  F.13,  if  the  total  number 
of  nodes  is  n. 

[k]  {q}  =  {f}  F.27 
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where 


T 

[kl    „    =    [L]   [C]  [L]  dV  F.28 

3nx3n 


tq}3nxl   =  tqxl  qy!  qzl  *  *  '  qxn  qyn  «W 

F.29 


m  T  T 

{f}  =       I  {P}        [N]    +  {t}        [N]     dA 

3nxl        i=j  i 

F.30 


In  this  study,  the  interpolation  functions  of  linear 
and  parabolic  isoparametric  elements  (eight  nodes  and 
twenty  nodes,  respectively)  [55]  are  chosen  for  shape 
functions.   Each  element  has  the  same  thickness  as  the  ply 
thickness  of  fiber  composite.   Once  the  nodal  displacement 
{q}  is  determined  by  solving  the  system  equations 
(Equation  F.27) ,  the  displacement  and  the  strain  fields 
can  be  obtained  by  Equations  F.23  and  F.24,  respectively. 
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